(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(id 



EP1 415 770 A1 



(12) 



EUROPEAN PATENT APPLICATION 

published in accordance with Art. 158(3) EPC 



(43) 


Date of publication: 


(51) 


Intel 7: B25J 5/00, B25J 13/00 




06.05.2004 Bulletin 2004/19 










(86) 


International application number: 


(21) 


Application number: 02736183.1 




PCT/JP2002/006467 


(22) 


Date of filing: 27.06.2002 


(87) 


International publication number: 








WO 2003/002309 (09.01 .2003 Gazette 2003/02) 


(84) 


Designated Contracting States: 


• 


IKEUCHI, Yasushi 




AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 




c/o K. K. Honda Gijyusu Kenkyusho 




MC NL PT SE TR 




Wako-shi, Saitama 351-0193 (JP) 




Designated Extension States: 


• 


KATO, Hisashi, 




AL LT LV MK RO SI 




c/o K. K. Honda Gijyutsu Kenkyusho 








4-1 , Chuo 1-chome (JP) 


(30) 


Priority 27.06.2001 US 300815 P 








13.07.2001 JP 2001214174 


(74) Representative: Hague, Alison Jane 




15.02.2002 JP 2002039201 




Frank B. Dehn & Co., 








European Patent Attorneys, 


(71) 


Applicant: Honda Glken Kogyo Kabushikl Kaisha 




179 Queen Victoria Street 




Tokyo 107-8556 (JP) 




London EC4V 4EL (GB) 


(72) 


Inventors: 






• 


KAWAI, Masakazu, 








c/o K. K. Honda Gtjyusu Kenkyusho 








Wako-shi, Saitama 351-0193 (JP) 







(54) METHOD OF ESTIMATING FLOOR REACTIONS OF BIPEDAL WALKING BODY; AND METHOD 
OF ESTIMATING JOINT MOMENTS OF BIPEDAL WALKING BODY 
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(57) It is judged whether a motion state of leg bodies 
2 is a single stance state or a double stance state. When 
it is the single stance state, a floor reaction force F acting 
on the leg body 2 touching the ground is estimated 
based on the equation of motion for the center of gravity 
GO of a bipedal walking moving body 1 . When it is the 
double stance state, floor reaction forces Fr, Ff respec- 
tively acting on the leg bodies 2 are estimated based on 
the equation of motion forthe center of gravity GO of the 
bipedal walking moving body 1, and the expression of 
relation between components of the floor reaction forces 
Fr, Ff on the respective leg bodies 2 and positions of 
ankle portions 1 2 of the leg bodies 2 relative to the cent- 
er of gravity GO of the bipedal walking moving body 1 . 
Further, based on an inverse dynamics model using es- 
timated values of the floor reaction forces Fr. Ff : mo- 
ments acting on a knee joint 1 0 and a hip joint 8 of each 
leg body 2 are estimated. Thereby, floor reaction forces 
acting on leg bodies of a bipedal walking moving body 
such as a human being, and moments acting on joints 
of the leg bodies are accurately grasped in real time by 
a relatively simple technique. 
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Description 

TECHNICAL FIELD 

5 [0001 ] The present invention relates to a method of estimating floor reaction forces acting on leg bodies of a bipedal 
walking moving body such as a human being or a bipedal walking robot, and further relates lo a method of estimating 
moments acting on joints of the leg bodies of the bipedal walking moving body by the use of estimated values of the 
floor reaction forces. 

w BACKGROUND ART 

[0002] For example, when executing an operation control of a walking assist device assisting walking motions of a 
human being, or a control of moving motions of a bipedal walking robot, it becomes necessary to successively grasp 
floor reaction forces acting on leg bodies (specifically, forces acting on ground touching portions of the leg bodies from 
15 a floor) of the human being or the bipedal walking robot. By grasping the floor reaction forces, it becomes possible to 
grasp moments or the like acting on joints of the leg bodies of the bipedal walking moving body and, based on the 
grasped moments or the like, it becomes possible to determine desired assist forces of the walking assist device, or 
desired driving torques for the respective joints of the bipedal walking robot or the like. 

[0003] As a technique of grasping the foregoing floor reaction forces, there has been known one disclosed in, for 
20 example, Laid-open Unexamined Patent Publication No. 2000-249570. According to this technique, in view of the fact 
that a waveform of temporal variation in floor reaction force on each leg body periodically changes during steady walking 
of a bipedal walking moving body, the floor reaction force on each leg body is grasped as a composite value (linear 
combination) of a plurality of trigonometric functions having mutually different periods being 1/n (n=1 , 2, *•■) of a walking 
period. In this case, as a weighting coefficient of each trigonometric function upon combining the plurality of trigono- 
us metric functions, there has been used a prescribed value determined in advance per bipedal walking moving body, or 
a value obtained by adjusting it depending on topography. 

[0004] However, inasmuch as the foregoing technique aims to grasp the floor reaction forces on the leg bodies with 
respect to one step or a plurality of steps of the bipedal walking moving body, when a gait of the bipedal walking moving 
body changes in sequence, it is difficult to accurately grasp the floor reaction forces. Further, for enhancing the accuracy 
30 of the floor reaction forces to be grasped, it is necessary to set the foregoing weighting coefficients of the trigonometric 
functions per bipedal walking moving body, or adjust them depending on topography or the like. Therefore, it is difficult 
to accurately grasp the floor reaction forces with a reduced influence of environment for movement of the bipedal 
walking moving body or a difference in individuality of the bipedal walking moving body. 

[0005] In case of a bipedal walking robot for example, there has been known one wherein force sensors such as 
35 six-axis force sensors are attached to ankle portions or foot portions of respective leg bodies, and floor reaction forces 

are grasped by outputs of these force sensors. Further, there has also been known a technique of walking a bipedal 

walking moving body on a force plate placed on a floor, and grasping floor reaction forces by outputs of the force plate. 

[0006] However, the technique using the force sensors has a disadvantage that when grasping floor reaction forces 

on leg bodies of particularly a human being, inasmuch as the force sensors should be attached to ankle portions or 
40 foot portions of the human being, the force sensors become obstructive to the walking in an ordinary living environment. 

Further, in case of the one using the force plate, the floor reaction forces can only be grasped in an environment where 

the force plate is arranged. 

[0007] The present invention has been made in view of the foregoing background and has an object to provide a 
floor reaction force estimating method that can accurately grasp floor reaction forces in realtime with a relatively simple 
45 technique, and that is suitable for grasping floor reaction forces with respect to particularly a human being as a bipedal 
walking moving body. 

[0008] Further, it is an object thereof to provide a method of estimating joint moments of the bipedal walking moving 
body, which can accurately grasp in real time the moments acting on joints such as knee joints of leg bodies by the 
use of estimated values of the floor reaction forces. 

50 

DISCLOSURE OF THE INVENTION 

[0009] At the outset, explanation will be given about a basic way of thinking about a floor reaction force estimating 
method of a bipedal walking moving body according to the present invention. 
55 [0010] As motion states of leg bodies of a bipedal walking moving body, for example, as motion states of the leg 
bodies upon walking, there are a single stance state wherein, as illustrated in Fig. 1 (a), only one leg body 2 (a leg body 
on a front side in an advancing direction in the figure) of both leg bodies 2, 2 of a bipedal walking moving body 1 touches 
the ground, and a double stance state wherein both leg bodies 2, 2 touch the ground as shown in Fig. 1 (b). 



3 

JNSDOCJD: <EP 1415770A1 I _> 



EP1 415 770 A1 



[0011] Here, first, in the foregoing single stance state, an equation of motion for the center of gravity (specifically, 
an equation of motion about translation of the center of gravity) of the bipedal walking moving body in an absolute 
coordinate system fixed relative to the ground on which the bipedal walking moving body moves is given by an expres- 
sion of relation such that the product of an acceleration of the center of gravity and a weight of the bipedal walking 

5 moving body is equal to a resultant force of gravity acting on the center of gravity (-the weight of the bipedal walking 
moving body x the acceleration of gravity) and a floor reaction force acting from the floor on a ground touching portion 
of the leg body touching the ground. Specifically, as shown in Fig. 1 (a), for example, in an absolute coordinate sysLem 
Cf fixed relative to a floor A, given that components of an acceleration a of the center of gravity GO of the bipedal 
walking moving body 1 in an X-axis direction (horizontal direction in the advancing direction of the bipedal walking 

10 moving body 1 ) and a Z-axis direction (vertical direction) are ax and az, respectively, and components of a floor reaction 
force F on the leg body 2 (leg body 2 on the supporting leg side) touching the ground in the X-axis direction and the 
Z-axis direction are Fx and Fz, respectively, the equation of motion for the center of gravity GO is expressed by the 
following equation (1 ). 

T (Fx, Fz-M * g) = M • T (ax, az) (1) 



(wherein M: weight of bipedal walking moving body, and g: acceleration of gravity) 
[0012] Parenthesis portions T ( , ) on both sides in the equation (1) represent two-component vectors. In the present 
specification, the notation in the form of T ( , ) represents a vector. 

[0013] Therefore, if the acceleration a= T (ax ; az) of the center of gravity GO of the bipedal walking moving body 1 is 
grasped, an estimated value of the floor reaction force F= T (Fx, Fz) can be derived by the following equation (2) using 
the acceleration a, a value of the weight M of the bipedal walking moving body 1 , and a value of the acceleration of 
gravity g. 

T (Fx, Fz) = M • T (ax, az-g) (2) 



[0014] In this case, the weight M necessary for deriving the estimated value of the floor reaction force F can be 
grasped in advance by measurement or the like. Further, with respect to a position of the center of gravity GO and the 
acceleration a, it is possible to grasp them in sequence according to a known technique or the like using outputs of 
sensors such as sensors for detecting bending angles (rotation angles) of the respective joints of the bipedal walking 
moving body 1 , accelero meters, gyro sensors, and so forth, which will be described later in detail. 
[0015] Next, an equation of motion for the center of gravity (specifically, an equation of motion about translation of 
the center of gravity) of the bipedal walking moving body in the foregoing both-leg ground touching state is given by 
an expression of relation such that the product of an acceleration of the center of gravity and the weight of the bipedal 
walking moving body is equal to a resultant force of gravity acting on the center of gravity (=the weight of the bipedal 
walking moving body x the acceleration of gravity) and floor reaction forces acting from the floor on the respective 
ground touching portions of both leg bodies (two floor reaction forces corresponding to both leg bodies, respectively). 
Specifically, as shown in Fig. 1(b), given that X and Z coordinate components of a floor reaction force Ff on the teg 
body 2 on the front side in the advancing direction of the bipedal walking moving body 1 are Ffx and Ffz, and X and Z 
coordinate components of a floor reaction force Fr on the leg body 2 on the rear side are Frx and Frz, the equation of 
motion for the center of gravity GO is expressed by the following equation (3). 

T (Ffx+Frx, Ffz+Frz-M • g) = M - T (ax, az) (3) 



[0016] The meaning of ax, az, M, and g in the equation (3) is as described before. 

[0017] On the other hand, according to knowledge of the inventors of the present application, in the double stance 
state, the floor reaction forces Ff, Fr respectively relating to the leg bodies 2, 2 can be deemed to almost act from 
specific portions 12f, 12r (e.g. ankle portions) in the vicinity of lower end portions of the respective leg bodies 2, 2 
toward the center of gravity GO of the bipedal walking moving body 1 , as shown in Fig. 1 (b). Then, there is established 
an expression of certain relation between a position of the specific portion 12f, 12r of the leg body 2 relative to the 
center of gravity GO, and the floor reaction force Ff, Fr acting on the leg body 2, that is, an expression of relation 
representing a relation wherein a direction of a segment Connecting between the center of gravity GO and the specific 
portion 12f, 12r of the leg body 2 (a direction of a position vector of the specific portion 12f, 12r relative to the center 

of gravity GO) is the same as a direction of the floor reaction force Ff, Fr relating to the leg body 2. 
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[0018] Specifically, referring to Fig. 1 (b) f given that coordinates of a position of the center of gravity GO in the foregoing 
absolute coordinate system Cf are (Xg, Zg), coordinates of a position of the specific portion 12f of the front-side leg 
body 2 are (Xf. Zf), and coordinates of a position of the specific portion 1 2r of the rear-side leg body 2 are (Xr. Zr), the 
foregoing expressions of relation become the following equations (4). 

5 

(Zf-Zg)/(Xf-Xg) - Ffz/Ffx 

(Zr-Zg)/(Xr-Xg) - Frz/Frx (4) 

10 

[0019] Then, the following equations (5) are obtained from these equations (4) and the foregoing equation (3). 

Ffx = M • {AXf . (AZr • ax - AXr • az 

- AXr • g)} /(AXf • AZr - AXr • AZf) 
Ffz ~ M • {AZf • (AZr • ax - AXr - az 
2Q - AXr • g)} /(AXf • AZr - AXr • AZf) 

Frx = M • {AXr • (-AZf - ax + AXf • az 

+ AXf - g)} /(AXf • AZr - AXr • AZf) 
Frz - M • {AZr - (-AZf - ax + AXf - az 

+ AXf • g)} /(AXf . AZr - AXr . AZf) (5) 

30 (wherein AZf=Xf-Xg, AZf=Zf-Zg, AXr=Xr-Xg, and AZr=Zr-Zg) 

[0020] Therefore, if an acceleration a= T (ax, az) of the center of gravity GO of the bipedal walking moving body 1 is 
grasped, and positions (these are expressed by AXf, AZf, AXr, and AZr in the equations (5)) of the specific portions 
12f, 12r respectively of the leg bodies 2, 2 relative to the center of gravity GO of the bipedal walking moving body 1 are 
grasped, estimated values of the floor reaction forces Ff= T (Ffx, Ffz), Fr= T (Frx, Frz) on the respective leg bodies 2, 2 

35 can be obtained by the foregoing equations (5) using a value of the weight M and a value of the acceleration of gravity 
g of the bipedal walking moving body 1 and the grasped acceleration a and the positions of the specific positions 1 2f , 1 2r. 
[0021] In this case, the weight M necessary for obtaining the estimated values of the floor reaction forces Ff, Fr can 
be grasped in advance by measurement or the like. With respect to the acceleration a of the center of gravity GO, the 
position of the center of gravity GO, and the positions of the specific portions 12f, 12r relative to the center of gravity 

40 GO, although details will be described later, it is possible to grasp them in sequence by a known technique or the like 
using outputs of the sensors detecting bending angles (rotation angles) of the respective joints of the bipedal walking 
moving body 1 , the accelerometers, the gyro sensors, and so forth. 

[0022] The present invention will be described next on the basis of what has bee n explained above. For accomplishing 
the foregoing objects, a method of estimating a floor reaction force for a bipedal walking moving body according to the 

45 present invention, i.e. a method of estimating a floor reaction force acting on each leg body of a bipedal walking moving 
body, comprises a first step of judging whether a motion state of the leg bodies of said bipedal walking moving body 
is a single stance state in which only one of the leg bodies touches the ground, or a double stance state in which both 
leg bodies touch the ground; a second step of sequentially deriving positions of the center of gravity of said bipedal 
walking moving body, and sequentially deriving accelerations of said center of gravity in an absolute coordinate system 

so fixed relative to the ground by the use of time series data about the positions of said center of gravity; and a third step 
of sequentially deriving positions of a specific portion relative to said center of gravity at least in said double stance 
state, said specific portion predetermined in the neighborhood of a lower end portion of each leg body. Then, the method 
of estimating a floor reaction force of the present invention is characterized by further comprising a step of deriving 
estimated values of said floor reaction force acting on the leg body touching the ground sequentially, in the single 

55 stance state of said bipedal walking moving body, based on an equation of motion for said center of gravity expressed 
by a weight and an acceleration of gravity of the bipedal walking moving body, the acceleration of said center of gravity, 
and said floor reaction force acting on the leg body touching the ground; and a step of deriving estimated values of 
said floor reaction forces respectively acting on both leg bodies sequentially, in the double stance state of said bipedal 
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walking moving body, based on an equation of motion for said center of gravity expressed by a weight and an accel- 
eration of gravity of the bipedal walking moving body, the acceleration of said center of gravity, and said floor reaction 
forces respectively acting on both leg bodies, and an expression of relation between the position of said specific portion 
of each leg body relative to said center of gravity and said floor reaction force acting on said leg body, said expression 
5 of relation being determined based on assumption that said floor reaction force acting on each leg body acts from said 
specific portion of said leg body toward said center of gravity. 

[0023] According to the present invention as described above, in the foregoing first step, it is judged whether the 
motion state of the leg bodies of the bipedal walking moving body is the single stance state or the double stance state, 
and the estimated values of the foregoing floor reaction forces are derived by the techniques depending on the respec- 
w tive supporting states. Specifically, in the single stance state of the bipedal walking moving body, the estimated value 
of the foregoing floor reaction force acting on the leg body touching the ground is derived from values of the weight, 
the acceleration of gravity, and the acceleration of the center of gravity of the bipedal walking moving body (see the 
equation (2)) based on the foregoing equation of motion for the center of gravity of the bipedal walking moving body 
(see the equation (1 )). 

15 [0024] On the other hand, in the double stance state of the bipedal walking moving body, the estimated values of 
the floor reaction forces acting on both leg bodies, respectively are derived from the weight, the acceleration of gravity, 
and the acceleration of the center of gravity of the bipedal walking moving body, and the position of the specific portion 
of each leg body relative to the center of gravity (see the equations (5)) based on the foregoing equation of motion for 
the center of gravity of the bipedal walking moving body (see the equation (3)) and the expression of relation (the 

20 equation (4)) between the position of the specific portion of each leg body relative to the foregoing center of gravity, 
and the foregoing floor reaction force acting on the subject leg body. Incidentally, in the single stance state, the floor 
reaction force acting on the free-leg side leg body (the leg body not touching the ground) is "0". 
[0025] In this case, the weight of the bipedal walking moving body necessary for deriving the estimated value of the 
floor reaction force may be grasped in advance by measurement or the like. With respect to the position of the center 

25 of gravity and the acceleration of the bipedal walking moving body, and the position of the specific portion of each leg 
body relative to the center of gravity, it is possible to grasp them in real time using data about outputs of sensors that 
are relatively small in size and easy to be attached to the bipedal walking moving body, such as sensors (potentiometers 
or the like) detecting bending angles (rotation angles) of the respective joints of the bipedal walking moving body, 
accelerometers, gyro sensors, and so forth. 

30 [0026] Therefore, according to the method of estimating a floor reaction force of the present invention, it is possible 
to grasp the floor reaction forces in real time by the relatively simple technique without attaching force sensors to ankle 
portions or foot portions of the bipedal walking moving body, or using a force place. 

[0027] In the method of estimating a floor reaction force of the present invention as described above, it is preferable 
that the specific portion of each leg body is an ankle portion of the leg body. Thereby, particularly, reliability of the 
35 foregoing assumption in the both-ieg ground touching state is enhanced. Therefore, the accuracy of not only the esti- 
mated value of the floor reaction force in the foregoing single stance state, but also the estimated values of the floor 
reaction forces in the double stance state can be enhanced. That is, the floor reaction force can be accurately estimated 
irrespective of the motion state of the leg bodies. 

[0028] Further, the method of estimating a floor reaction force of the present invention comprises a step of measuring 
40 an acceleration, in an upward/downward direction, of a lower portion of a body supported on both leg bodies via a hip 
joint of each leg body, said lower portion of the body located near said hip joints, wherein, in said first step, the motion 
state of said bipedal walking moving body is judged such that when the acceleration of said lower portion of the body 
in the upward/downward direction increases to a predetermined threshold value or more, said double stance state 
starts while said single stance state finishes and, when the estimated value of said floor reaction force acting on the 
45 leg body which precedingly makes a takeoff is lowered to a predetermined threshold value or less in said double stance 
state, said double stance state finishes while said single stance state starts. 

[0029] Specifically, upon transition of the motion state of the leg bodies from the single stance state to the double 
stance state during movement (walking) of the bipedal walking moving body, landing of the free-leg side leg body 
causes the acceleration of a lower part of a body in an upward/downward direction (upward acceleration) to temporarily 

50 become remarkably large. This phenomenon is not generated normally in the other motion state of the leg bodies. On 
the other hand, upon transition of the motion state of the leg bodies from the double stance state to the single stance 
state, a takeoff operation of one of the leg bodies causes the floor reaction force acting on that leg body to be lowered 
to "0". Therefore, by judging the motion state of the leg bodies as described above, judgement as to whether it is the 
single stance state or the double stance state can be properly carried out. As a result, a floor reaction force estimated 

55 value calculation method that differs between the single stance state and the double stance state can be switched at 
appropriate timing to thereby increase the accuracy of the estimated value of the floor reaction force. The acceleration, 
in the upward/downward direction, of the lower part of the body necessary for judging the motion state of the leg bodies 
can be easily grasped from an output of an accelerometer by, for example, attaching the accelerometer to the lower 
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part of the body. 

[0030] Like a human being, if the foregoing body has a waist coupled to both leg bodies via hip joints, and a chest 
located on the waist so as to be tilted freely relative to the waist it is preferable that the acceleration, in the upward/ 
downward direction, of the lower part of the body to be measured is an acceleration of the waist in the upward/downward 
5 direction. 

[0031] Further, in Lhe method o\ estimating a floor reaction force of the present invention, various techniques are 
considered as a technique of grasping the position of the center of gravity of the bipedal walking moving body and the 
acceleration of the center of gravity in the foregoing second step, and it is possible to use various known methods. 
However, it is preferable to grasp the position of the center of gravity and the acceleration of the center of gravity 

10 according to the following method. 

[0032] Specifically, the method comprises a step of respectively measuring an inclination angle of a body supported 
on both leg bodies via a hip joint of each leg body, bending angles respectively of at least the hip joint and a knee joint 
of each leg body, and an acceleration of a predetermined reference point of said bipedal walking moving body in said 
absolute coordinate system, wherein, in said second step, based on the inclination angle of said body, the bending 

15 angles respectively of said hip joints and said knee joints, a rigid body link model formed by expressing said bipedal 
walking moving body as a linked body of a plurality of rigid bodies, prederived weights of respective rigid-body corre- 
sponding portions of the bipedal walking moving body corresponding to the respective rigid bodies of said rigid body 
link model, and positions of the prederived centers of gravity of the rigid-body corresponding portions in the respective 
rigid-body corresponding portions, positions of the center of gravity of said bipedal walking moving body relative to 

20 said reference point are sequentially derived, accelerations of said center of gravity relative to said reference point are 
sequentially derived based on time series data about the positions of said center of gravity, and the acceleration of 
said center of gravity in said absolute coordinate system is derived from the acceleration of said center of gravity 
relative to said reference point, and the acceleration of said reference point in said absolute coordinate system. 
[0033] Specifically, when the reference point is arbitrarily set to the bipedal walking moving body, the position of the 

25 center of gravity of the bipedal walking moving body relative to the reference point is almost determined by a mutual 
posture relationship among the body, a thigh portion of each leg body from the hip joint to the knee joint, and a cms 
portion on a lower side below the knee joint. By measuring the inclination angle of the body and the bending angles 
respectively of the hip joints and the knee joints, the posture relationship can be grasped from those measurement 
data. Further, although details will be described later, when the foregoing rigid body link model is supposed for the 

30 bipedal walking moving body (e.g. a model regarding as rigid bodies a portion (including the body) on an upper side 
of the hip joints of both leg bodies of the bipedal walking moving body, and the thigh portion and the crus portion of 
each leg body) , the position of the center of gravity of the bipedal walking moving body relative to thef oregoing reference 
point can be derived based on the weights of the respective rigid-body corresponding portions of the bipedal walking 
moving body, the positions of the centers of gravity of the rigid-body corresponding portions in the respective rigid- 

35 body corresponding portions (specifically, the positions of the rigid-body corresponding portions in coordinate systems 
fixed to the respective rigid-body corresponding portions), and the foregoing posture relationship. Further, the accel- 
eration of the center of gravity relative to the reference point can be derived as a second-order differentiated value of 
the position of the center of gravity which is grasped from time series data about the positions of the center of gravity. 
Therefore, by measuring the acceleration of the reference point in the foregoing absolute coordinate system, the ac- 

40 celeration of the center of gravity of the bipedal walking moving body in the absolute coordinate system can be derived 
as a composite acceleration of the acceleration of the center of gravity relative to the reference point and the acceler- 
ation of the reference point. 

[0034] In this case, the inclination angles of the body necessary for grasping the acceleration of the bipedal walking 
moving body as described above can be grasped from outputs of accelero meters and gyro sensors, or inclinometers 

45 or the like attached to the body, and the bending angles of the hip joint and the knee joint of each leg body can be 
grasped from outputs of sensors such as potentiometers attached to the portions of the respective joints. Further, the 
acceleration of the reference point in the absolute coordinate system can be grasped from outputs of sensors such as 
accelerometers attached to a portion integral with the reference point. Further, the weights of the respective rigid-body 
corresponding portions of the bipedal walking moving body and the positions of the centers of gravity of the rigid-body 

50 corresponding portions in the respective rigid-body corresponding portions can be grasped in advance by measurement 
or the like. 

[0035] Therefore, it becomes possible to easily grasp in real time the position and acceleration of the center of gravity 
of the bipedal walking moving body without mounting the relatively large sensors and so forth onto the bipedal walking 
moving body. 

55 [0036] When deriving the position and acceleration of the center of gravity of the bipedal walking moving body as 
described above, it is preferable that the reference point is set to the body. Thereby, the sensors such as the acceler- 
ometers for measuring the acceleration of the reference point in the absolute coordinate system can be mounted on 

the body, and therefore, the sensors attached to the leg bodies can be reduced, thereby to maKe it possible to avoid 
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the situation where those sensors become obstructive to walking motions of the bipedal walking moving body. 
[0037] On the other hand, like a human being, if the body has a waist coupled to both leg bodies via hip joints, and 
a chest located on the waist so as to be tiltable relative to the waist, it is preferable that the inclination angle of the 
body used for deriving the position of the center of gravity comprises inclination angles respectively of the waist and 
5 the chest. Particularly in this case, it is preferable that the foregoing rigid body link model is a model expressing a crus 
portion located on a lower side of the knee joint of each leg body of the bipedal walking moving body, a thigh portion 
between the knee joint and the hip joint, the waist, and an upper body portion located on an upper side of the waist 
and including the chest, as rigid bodies, respectively. 

[0038] In accordance therewith, particularly when the bipedal walking moving body is a human being, the position 
10 of the center of gravity thereof and the acceleration thereof can be grasped accurately to thereby enhance the accuracy 
of the estimated values of the floor reaction forces. 

[0039] When deriving the position of the center of gravity of the bipedal walking moving body relative to the reference 
point based on the inclination angles of the body and so forth as described above, the position of the foregoing specific 
portion of each leg body relative to the center of gravity of the bipedal walking moving body can be derived by deriving 

15 the position of the specific portion relative to the reference point in the foregoing third step. In this case, parameters 
necessary for deriving the position of the specific portion relative to the reference point differ depending on a set position 
of the reference point. For example, when the specific portion is an ankle portion of each leg body, and the reference 
point is set to the body, it is possible to derive the position of the specific portion of each leg body relative to the reference 
point based on the inclination angle of the body (the inclination angle of the waist when the body has the waist and the 

20 chest, and the reference point is set to the waist), the bending angles of the hip joint and the knee joint of each leg 
body, and the sizes (lengths) of the thigh portion and the crus portion of each leg body. 

[0040] Next, a method of estimating a moment for a bipedal walking moving body according to the present invention 
is a method of estimating a moment acting on at least one joint of each leg body of the bipedal walking moving body 
by the use of the estimated values of the floor reaction force on each leg body sequentially derived by the foregoing 

25 method of estimating a floor reaction force. Then, the method of estimating a moment of the present invention is- 
characterized by comprising a step of respectively measuring an inclination angle of a body supported on both leg 
bodies via a hip joint of each leg body, bending angles respectively of at least the hip joint and a knee joint of each leg 
body, and an acceleration of a predetermined reference point of said bipedal walking moving body in said absolute 
coordinate system; a step of sequentially deriving inclination angles of respective rigid-body corresponding portions 

30 of the bipedal walking moving body corresponding to the respective rigid bodies of said rigid body link model, based 
on the inclination angle of said body, the bending angles respectively of said hip joint and said knee joint of each leg 
body, and a rigid body link model formed by expressing said bipedal walking moving body as a linked body of a plurality 
of rigid bodies; a step of. sequentially deriving positions of the centers of gravity of the respective rigid-body corre- 
sponding portions relative to said reference point, based on the inclination angles of said respective rigid-body corre- 

35 sponding portions, prederived weights of the respective rigid-body corresponding portions, and positions of the pred- 
erived centers of gravity of the rigid-body corresponding portions in the respective rigid-body corresponding portions, 
and sequentially deriving accelerations of the centers of gravity of the respective rigid-body corresponding portions 
relative to said reference point based on time series data about the positions of the centers of gravity of the respective 
rigid-body corresponding portions; a step of sequentially deriving accelerations of the centers of gravity of the respective 

40 rigid-body corresponding portions in said absolute coordinate system from the accelerations of the centers of gravity 
of the respective rigid-body corresponding portions relative to said reference point, and the accelerations of said ref- 
erence point in said absolute coordinate system; a step of sequentially deriving angular velocities of the respective 
rigid-body corresponding portions based on time series data about the inclination angles of the respective rigid-body 
corresponding portions; a step of sequentially deriving estimated positions of a floor reaction force acting point of each 

45 leg body in said bipedal walking moving body, based on at least one of an inclination angle of a thigh portion of said 
leg body and the bending angle of the knee joint of said leg body as the rigid-body corresponding portions of said 
bipedal walking moving body; and a step of estimating a moment acting on at least one of the joints of each leg body 
of said bipedal walking moving body based on an inverse dynamics model using the estimated value of said floor 
reaction force, the estimated position of said floor reaction force acting point, the accelerations of the centers of gravity 

50 of the respective rigid-body corresponding portions and the angular velocities of said rigid-body corresponding portions 
in said absolute coordinate system, the inclination angles of the respective rigid-body corresponding portions, the 
prederived weights and sizes of the respective rigid-body corresponding portions, the positions of the prederived cent- 
ers of gravity of the rigid-body corresponding portions in the respective rigid-body corresponding portions, and pred- 
erived moments of inertia of the respective rigid-body corresponding portions. 

55 [0041 ] When, as described before, the second step in the foregoing floor reaction force estimating method comprises 
a step of respectively measuring the inclination angle of the body ; the bending angles respectively of the hip joint and 
the knee joint of each leg body, and the acceleration of the reference point of the bipedal walking moving body in the 
absolute coordinate system, for deriving the position of the center of gravity of the bipedal walking moving body relative 
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to the reference point, and so forth, it is not necessary to newly measure them. Further the rigid body link model may 
be the same as the rigid body link model used for deriving the position of the center of gravity of the bipedal walking 
moving body, and so forth. 

[0042] In the method of estimating a moment of the present invention as described above, by measuring the incli- 
nation angle of the body and the bending angles respectively of the hip joints and the knee joints, it is possible to grasp, 
from those measurement data, inclination angles of the respective rigid-body corresponding portions (these represent 
a mutual posture relationship of the respective rigid-body corresponding portions), such as the body, the thigh portions, 
the cms portions, and so forth, of the bipedal walking moving body. Then, based on the weights of the respective rigid- 
body corresponding portions, the positions of the centers of gravity of the rigid-body corresponding portions in the 
respective rigid-body corresponding portions (specifically, the positions of the rigid-body corresponding portions in 
coordinate systems fixed to the respective rigid-body corresponding portions), and the inclination angles of the respec- 
tive rigid-body corresponding portions, it is possible to derive the positions of the centers of gravity of the respective 
rigid-body corresponding portions with respect to the reference point. Further, the acceleration of the center of gravity 
of each rigid-body corresponding portion relative to the reference point is derived as a second-order differentiated 
value of the position of the center of gravity which is grasped from time series data about the positions of the center 
of gravity. Therefore, by measuring the acceleration of the reference point in the absolute coordinate system, the 
acceleration of the center of gravity of each rigid-body corresponding portion of the bipedal walking moving body in 
the absolute coordinate system can be derived as a composite acceleration of the acceleration of the center of gravity 
relative to the reference point and the acceleration of the reference point (the acceleration in the absolute coordinate 
system). 

[0043] Further, an angular velocity of each rigid-body corresponding portion is derived as a second-order differenti- 
ated value of the inclination angle which is grasped from time series data about the inclination angles of each rigid- 
body corresponding portion. 

[0044] According to knowledge of the inventors of the present application, the position of a floor reaction force acting 
point of each leg body in the bipedal walking moving body, for example, the position of a floor reaction force acting 
point of each leg body relative to the ankle portion of the leg body, has a close correlation with the inclination angle of 
the thigh portion of each leg body as the rigid-body corresponding portion of the bipedal walking moving body, and the 
bending angle of the knee joint of each leg body. Therefore, it is possible to derive an estimated position of the floor 
reaction force acting point in the bipedal walking moving body based on at least one of the inclination angle of the thigh 
portion of each leg body and the bending angle of the knee joint thereof. 

[0045] Then, when the accelerations of the centers of gravity of the respective rigid-body corresponding portions of 
the bipedal walking moving body, the angular velocities of the respective rigid-body corresponding portions, and the 
estimated position of the floor reaction force acting point are derived as described above, it is possible to estimate 
moments acting on the knee joint and the hip joint of each leg body based on a known so-called inverse dynamics 
model, using those data and so forth along with the estimated values of the floor reaction forces derived by the foregoing 
floor reaction force estimating method. Briefly speaking, the technique based on this inverse dynamics model uses the 
equation of motion for the translation motion of the center of gravity of each rigid-body corresponding portion of the 
bipedal walking moving body, and the equation of motion forthe rotating motion of the rigid-body corresponding portion 
(e.g. the rotating motion about the center of gravity of the rigid-body corresponding portion), and derives moments 
acting on the respective joints of the bipedal walking moving body corresponding to the respective joints of the rigid 
body link model, in order from the joint closer to the floor reaction force acting point. Although details will be described 
later, assuming that, for example, each leg body is a linked body having a thigh portion and a crus portion as rigid- 
body corresponding portions, respectively, a force (joint reaction force) acting on a knee joint of the leg body can be 
grasped by applying an acceleration of the center of gravity of the crus portion, an estimated value of a floor reaction 
force acting on the leg body, and a value of the weight of the crus portion, into the equation of motion for the translation 
motion of the center of gravity of the crus portion of each leg body. Further, a moment on the knee joint of the leg body 
can be estimated by applying a joint reaction force acting on the knee joint of the leg body, an angular velocity of the 
crus portion of the leg body, an estimated position of a floor reaction force acting point of the leg body, an estimated 
value of a floor reaction force on the leg body, data values relating to a position of the center of gravity of the crus 
portion in the crus portion and a size (length) of the crus portion, a value of a moment of inertia of the crus portion, and 
a value of an inclination angle of the crus portion, into the equation of motion forthe rotating motion of the crus portion. 
[0046] Further, a joint reaction force acting on a hip joint of the leg body can be grasped by applying an acceleration 
of the center of gravity of the thigh portion, a joint reaction force acting on the knee joint of the leg body, and a value 
of the weight of the thigh portion, into the equation of motion for the translation motion of the center of gravity of the 
thigh portion. Further, a moment on the hip joint of the leg body can be estimated by applying joint reaction forces 
respectively acting on the knee joint and the hip joint of the leg body, an angular velocity of the thigh portion of the leg 
body, data values relating to a position of the center of gravity of the thigh portion in the thigh portion and a size (length) 
of the thigh portion, a value of a moment of inertia of the thigh portion, and a value of an inclination angle of the thigh 
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portion, into the equation of motion for the rotating motion of the thigh portion. 

[0047] According to the method of estimating a moment of the present invention as described above, by estimating 
the moments acting on the joints of the leg bodies using the estimated values of the floor reaction forces derived by 
the foregoing floor reaction force estimating method of the present invention, it is possible to relatively accurately 
5 estimate in real time the moments acting on the joints of the leg bodies without attaching the relatively large sensors 
and so forth to the bipedal walking moving body. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 [0048] Fig. 1 is a diagram for explaining the basic principle of a floor reaction force estimating method of the present 
invention, Fig. 2 is a diagram exemplarily showing a human being as a bipedal walking moving body and structures of 
devices attached to the human being in one embodiment of the present invention, Fig. 3 is a block diagram for explaining 
functions of an arithmetic processing unit provided in the devices of Fig. 2. and Fig. 4 is a diagram showing a rigid 
body link model used in the processing of the arithmetic processing unit of Fig. 3. Fig. 5 is a diagram showing a 

15 correlation between advancing direction components of a floor reaction force acting point vector and an inclination 
angle of a thigh portion during normal walking, Fig. 6 is a diagram showing a correlation between vertical direction 
components of a floor reaction force acting point vector and an inclination angle of a thigh portion during normal walking, 
Fig. 7 is a diagram showing a correlation between advancing direction components of a floor reaction force acting point 
vector and a bending angle of a knee joint upon sitting in a chair, Fig. 8 is a diagram showing a correlation between 

20 advancing direction components of a floor reaction force acting point vector and a bending angle of a knee joint upon 
standing up from a chair, Fig. 9 is a diagram showing a correlation between advancing direction components of a floor 
reaction force acting point vector and an inclination angle of a thigh portion upon going upstairs, and Fig. 1 0 is a diagram 
showing a correlation between advancing direction components of a floor reaction force acting point vector and an 
inclination angle of a thigh portion upon going downstairs. Fig. 11 is a diagram for explaining the processing in joint 

25 moment estimating means of the arithmetic processing unit of Fig. 3, Fig. 12 is a graph illustrating the states of temporal 
variations of estimated values of floor reaction forces during normal walking, which are derived by the embodiment of 
the present invention, Fig. 1 3 is a graph illustrating the states of temporal variations of estimated values of floor reaction 
forces upon going upstairs, which are derived by the embodiment of the present invention, Fig. 1 4 is a graph illustrating 
the state of temporal variation of estimated values of floor reaction forces upon standing up from a chair, which are 

30 derived by the embodiment of the present invention, and Fig. 15 is a graph illustrating the states of temporal variations 
of estimated values of moments on a knee joint and a hip joint, which are derived by the embodiment of the present 
invention. 

BEST MODE FOR CARRYING OUT THE INVENTION 

35 

[0049] One embodiment of the present invention will be described with reference to the foregoing figure 1 and Figs. 
2 to 1 0. 

[0050] The present embodiment is an embodiment wherein a floor reaction force estimating method and a method 
of estimating a moment of the present invention are applied to a human being as a bipedal walking moving body. 

40 [0051] As exemplarily shown in Fig. 2. a human being 1, when its structure is roughly classified, has a pair of left 
and right leg bodies 2, 2, a body 5 comprising a waist 3 and a chest 4, a head 6, and a pair of left and right arm bodies 
7, 7. The waist 3 is coupled to the respective leg bodies 2, 2 via a pair of left and right hip joints 8, 8 so that the body 
5 is supported on both leg bodies 2, 2. Further the chest 4 of the body 5 is disposed on an upper side of the waist 3 
so as to be free inclined toward a forward side of the human being 1 relative to the waist 3. And, the arm bodies 7, 7 

45 extend from both left and right side portions of an upper part of the chest 4, and the head 6 is supported at an upper 
end portion of the chest 4. 

[0052] Each of the leg bodies 2, 2 has a thigh portion 9 extending from the hip joint 8, and a crus portion 1 1 extending 
from a tip of the thigh portion 9 via a knee joint 1 0, and a foot portion 1 3 is coupled to a tip portion of the crus portion 
11 via an ankle portion (ankle joint) 12. 
50 [0053] In the present embodiment, for estimating a floor reaction force acting on each leg body 2 of the human being 
1 having the foregoing structure, and further estimating moments acting on the knee joint 10 and the hip joint 8, the 
human being 1 is equipped with the following devices. 

[0054] Specifically, the chest 4 of the body 5 is mounted with a gyro sensor 14 (hereinafter referred to as the chest 
gyro sensor 14) that produces an output corresponding to an angular velocity following an inclination of the chest 4, 
55 an accelerometer 1 5 (hereinafter referred to as the chest forward/backward accelerometer 1 5) that produces an output 
corresponding to an acceleration of the chest 4 in a forward/backward direction, an arithmetic processing unit 16 com- 
prising a CPU, a RAM, a ROM, and so forth, and a battery 17 serving as a power supply for the arithmetic processing 
unit 16 and so forth. In this case, these chest gyro sensor 14, chest forward/backward accelerometer 15, arithmetic 
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processing unit 1 6, and battery 1 7 are : for example, received in a reception member 1 8 in the form of a shoulder bag 
that is fixed to the chest 4 using a belt or the like not shown, and unitariiy fixed to the chest 4 via the reception member 18. 
[0055] The acceleration represented by the output of the chest accelerometer 1 5 is : more specifically, an acceleration 
in a forward/backward direction in a direction of horizontal section of the chest 4 (a direction perpendicular to a central 

5 axis of the chest 4). In the state where the human being 1 is standing erect on a flatland. it is an acceleration in the 
forward/backward horizontal direction (X-axis direction of the absolute coordinate system in Fig. 2). while, in the stale 
where the waist 3 or the chest 4 is inclined from the vertical direction (Z-ax^s direction of the absolute coordinate system 
in Fig. 2), it is an acceleration in a direction inclined relative to the horizontal direction by an inclination angle of the 
chest 4 relative to the vertical direction. 

10 [0056] Further a gyro sensor 1 9 (hereinafter referred to as the waist gyro sensor 1 9) for producing an output corre- 
sponding to an angular velocity following an inclination of the waist 3, an accelerometer 20 (hereinafter referred to as 
the waist forward/backward accelerometer 20) for producing an output corresponding to an acceleration of the waist 
3 in a forward/backward direction, and an accelerometer 21 (hereinafter referred to as the waist upward/downward 
accelerometer 21) for producing an output corresponding to an acceleration of the waist 3 in an upward/downward 

15 direction are unitariiy mounted/fixed to the waist 3 of the body 5 via fixing means such as a belt not shown. 

[0057] Here, the waist forward/backward accelerometer 20 is, more specifically, a sensor, like the chest forward/ 
backward accelerometer 15, for detecting an acceleration in a forward/backward direction in a direction of horizontal 
section of the waist 3 (a direction perpendicular to a central axis of the waist 3). Further, the waist upward/downward 
accelerometer 21 is. more specifically, a sensor for detecting an acceleration in an upward/downward direction in a 

20 direction of the central axis of the waist 3 (this is perpendicular to the acceleration detected by the waist forward/ 
backward accelerometer 20). Incidentally, the waist forward/backward accelerometer 20 and the waist upward/down- 
ward accelerometer 21 may be unitariiy formed as a two-axis accelerometer. 

[0058] Further the hip joint 8 and the knee joint 1 0 of each leg body 2 are mounted with a hip joint angle sensor 22 
and a knee joint angle sensor 23 that produce outputs corresponding to respective bending angles Atic and AGd. With 
25 respect to the hip joint angle sensors 22, although only the hip joint angle sensor 22 relating to the hip joint 8 of the 
leg body 2 on this side (the right side of the human being 1 facing forward) is shown in Fig. 2 ; the hip joint 8 of the leg 
body 2 on the other side (the left side of the human being 1 facing forward) is mounted with the hip joint angle sensor 
22 concentrically with the hip joint angle sensor 22 on this side. 

[0059] These angle sensors 22 and 23 are formed by : for example, potentiometers, and attached to each leg body 

30 2 via means such as a band member not shown. Here, the bending angle A0c detected by each hip joint angle sensor 
22 is, more specifically, a rotation angle of the thigh portion 9 of the leg body 2 about the hip joint 8 (about a central 
axis of the hip joint 8 in a leftward/rightward direction of the human being 1 ) relative to the waist 3, using as a criterion 
the state where a posture relationship between the waist 3 and the thigh portion 9 of the leg body 2 is in a predetermined 
posture relationship (e.g. a posture relationship in which the central axis of the waist 3 and a central axis of the thigh 

35 portion 9 become substantially parallel to each other as in an erect posture state of the human being 1). Likewise, the 
bending angle A6d detected by each knee joint angle sensor 23 is a rotation angle of the crus portion 11 about the 
knee joint 10 (about a central axis of the knee joint 10 in a leftward/rightward direction of the human being 1) relative 
to the thigh portion 9, using as a criterion the state where a posture relationship between the thigh portion 9 and the 
crus portion 11 of the leg body 2 is in a predetermined posture relationship (e.g. a posture relationship in which the 

40 central axis of the thigh portion 9 and a central axis of the crus portion 11 become substantially parallel to each other). 
[0060] Incidentally, the foregoing respective sensors 14, 15, and 19 to 23 are connected to the arithmetic processing 
unit 16 via signal lines not shown, for inputting outputs thereof into the arithmetic processing unit 16. 
[0061] The foregoing arithmetic processing unit 1 6 is provided with functional means as shown in Fig. 3. Specifically, 
the arithmetic processing unit 16 is provided with leg body motion judging means 24 forjudging whether the motion 

45 state of the leg bodies 2, 2 of the human being 1 is the single stance state (state of Fig. 1(a)) or the double stance 
state (state of Fig. 1(b)), using detection data of the waist upward/downward accelerometer 21 and data about an 
estimated value of a floor reaction force of each leg body 2 derived by later-described floor reaction force estimating 
means 36. Further, the arithmetic processing unit 16 is provided with chest inclination angle measuring means 25 for 
measuring an inclination angle 6a of the chest 4 in the absolute coordinate system Cf (specifically, an inclination angle 

50 ea relative to, for example, the vertical direction; see Fig. 2), using detection data of the chest forward/backward ac- 
celerometer 1 5 and the chest gyro sensor 1 4, and waist inclination angle measuring means 26 for measuring an incli- 
nation angle 9b of the waist 3 in the absolute coordinate system Cf (specifically, an inclination angle 0b relative to, for 
example, the vertical direction; see Fig. 2), using detection data of the waist forward/backward accelerometer 20 and 
the waist gyro sensor 1 9. 

55 [0062] Further, the arithmetic processing unit 16 is provided with reference acceleration measuring means 27 for 
deriving an acceleration (translation acceleration) a 0 = T (a 0 x ; aoZ), in the absolute coordinate system Cf, of the origin O 
of a body coordinate system Cp (xz coordinates in Fig. 2) which is set at the waist 3 as a reference point of the human 
being 1 in the present embodiment as shown in Fig. 2, using detection data of the waist forward/backward accelerometer 
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20 and the waist upward/downward accelerometer 21 and data about the inclination angle 9b of the waist 3 measured 
by the foregoing waist inclination angle measuring means 26. Here, the body coordinate system Cp is, more specifically; 
a coordinate system (coordinate system of which directions of three axes are the same as those of the foregoing 
absolute coordinate system Cf) wherein the origin O is set to. for example, a middle point of a line connecting the 

5 respective centers of the left and right hip joints 8, 8 of the human being 1 , a z-axis direction is set to the vertical 
direction, and an x-axis direction is set to the horizontal direction headed forward of the human being 1 . 
[0063] Further, the arithmetic processing unit 1 6 is provided with leg body posture calculating means 28 for deriving 
inclination angles 6c and 9d respectively of the thigh portion 9 and the crus portion 1 1 of each leg body 2 in the absolute 
coordinate system Cf (specifically, inclination angles 6c and 9d relative to. for example, the vertical direction: see Fig. 

10 2). using detection data of the hip joint angle sensor 22 and the knee joint angle sensor 23 of the leg body 2 and data 
about the inclination angle Bb of the waist 3 by the foregoing waist inclination angle measuring means 26. 
[0064] Further, the arithmetic processing unit 1 6 is provided with respective-portion center of gravity position calcu- 
lating means 29 for deriving positions of the centers of gravity of respective rigid-body corresponding portions (specif- 
ically, positions of the centers of gravity of respective rigid-body corresponding portions in the foregoing body coordinate 

15 system Cp) of the human being 1 corresponding to a later-described rigid body link model, using data about the incli- 
nation angle 9a of the chest 4, the inclination angle 6b of the waist 3, and the inclination angle 9c of the thigh portion 
9 and the inclination angle Gd of the crus portion 1 1 of each leg body 2 obtained by the foregoing chest inclination angle 
measuring means 25, waist inclination angle measuring means 26, and leg body posture calculating means 28, body 
center of gravity position calculating means 30 for deriving a position of the center of gravity of the whole human being 

20 1 in the foregoing body coordinate system Cp, using data about the positions of the centers of gravity of the respective 
rigid-body corresponding portions, ankle position calculating means 31 for deriving, using data about the position of 
the center of gravity GO of the whole human being 1 (see Fig. 1 ; hereinafter referred to as the body center of gravity 
GO) and data about the respective inclination angles 9c and 9d of the thigh portion 9 and the crus portion 11 of each 
body leg 2 by the foregoing leg body posture calculating means 28, a position of the ankle portion 1 2 of each leg body 

25 2, as a specific part of each leg 2 in the present embodiment, relative to the body center of gravity GO (specifically, 
AXf, AZf, AXr, and AZr in the foregoing equations (5)), and body center of gravity acceleration calculating means 32 
for deriving an acceleration a= T (ax, az) (see Fig. 1) of the body center of gravity GO in the absolute coordinate system 
Cf . using data about the position of the body center of gravity by the foregoing body center of gravity position calculating 
means 30 and data about the acceleration a 0 of the origin O of the body coordinate system Cp by the foregoing reference 

30 acceleration measuring means 27. 

[0065] Further, the arithmetic processing unit 1 6 is provided with leg body respective-portion acceleration calculating 
means 33 for deriving accelerations (translation accelerations) of the centers of gravity of the thigh portion 9 and the 
crus portion 11 , respectively, of each leg body 2 in the absolute coordinate system Cf, using data about the positions 
of the centers of gravity of the respective rigid-body corresponding portions of the human being 1 (specifically, the 

35 positions of the centers of gravity of the rigid-body corresponding portions relating to each leg body 2) by the foregoing 
respective-portion center of gravity position calculating means 29 and data about the acceleration a 0 of the origin O 
of the body coordinate system Cp by the foregoing reference acceleration measuring means 27, leg body respective- 
portion angular velocity calculating means 34 for deriving angular velocities of the thigh portion 9 and the crus portion 
11 of each leg body 2 in the absolute coordinate system Cf, using data about the respective inclination angles 9c and 

40 ed of the thigh portion 9 and the crus portion 1 1 of the leg body 2 by the foregoing leg body posture calculating means 
28, and floor reaction force acting point estimating means 35 for estimating a position of a floor reaction force acting 
point of each leg body 2 touching the ground, using data about the inclination angle 9c of the thigh portion 9 derived 
by the foregoing leg body posture calculating means 28 and the bending angle A9d of the knee joint 10 measured by 
the foregoing knee joint angle sensor 23. 

45 [0066] Further, the arithmetic processing unit 1 6 is provided with floor reaction force estimating means 36 for deriving 
an estimated value of a floor reaction force acting on each leg body 2, using data about the acceleration a of the body 
center of gravity by the foregoing body center of gravity acceleration calculating means 32, data about the position of 
the ankle portion 1 2 of each leg body 2 relative to the body center of gravity by the foregoing ankle position calculating 
means 31 , and data about the judgement result on the motion state of the leg bodies 2, 2 by the foregoing leg body 

50 motion judging means 24, and joint moment estimating means 37 for estimating moments acting on the knee joint 10 
and the hip joint 8, respectively, of each leg body 2, using data about this estimated value of the floor reaction force, 
data about the accelerations of the centers of gravity of the thigh portion 9 and the crus portion 1 1 of each leg body 2 
by the foregoing leg body respective -portion acceleration calculating means 33, data about the angular velocities of 
the thigh portion 9 and the crus portion 11 of each leg body 2 by the leg body respective-portion angular velocity 

55 calculating means 34, data about the estimated position of the floor reaction force acting point by the floor reaction 
force acting point estimating means 35, and data about the inclination angles 9c and 9d of the thigh portion 9 and the 
crus portion 11 , respectively, of each leg body 2 by the foregoing leg body posture calculating means 28. 
[0067] Now, an operation of this embodiment will be described along with more detailed processing contents of the 
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respective means of the foregoing arithmetic processing unit 16. 

[0068] In the present embodiment, when the human being 1 turns on a power supply switch, not shown, of the 
arithmetic processing unit 1 6 in the state where both leg bodies 2, 2 are on the floor, upon, for example, moving the 
leg bodies 2, 2, such as walking, the processing by the arithmetic processing unit 16 is executed in sequence per 
5 predetermined cycle time as will be described hereinbelow, to thereby successively derive estimated values of floor 
reaction forces and so forth acting on each leg body 2. 

[0069] Specifically, at the outset, the arithmetic processing unit 1 6 executes the processing of the foregoing leg body 
motion judging means 24. In the processing of the leg body motion judging means 24, detection data about the accel- 
eration of the waist 3 in the upward direction by the foregoing waist upward/downward accelerometer 21 is compared 

10 with a prescribed threshold value determined in advance per the foregoing cycle time. Then, when a detected value 
of the acceleration exceeds the threshold value, it is judged that the double stance state as shown in the foregoing 
figure 1 (b) starts, while the single stance state as shown in the foregoing figure 1 (a) finishes. Specifically, upon transition 
from the single stance state to the double stance state during walking of the human being 1 , landing of the free-ieg 
side leg body 2 on the floor causes the waist 3 near the hip joints 8 to produce a substantially upward and relatively 

15 large acceleration (acceleration that can not be produced in the normal single stance state). Therefore, by comparing 
the detection data about the acceleration of the waist 3 in the upward direction by the waist upward/downward accel- 
erometer 21 with the prescribed threshold value as described above, the foregoing leg body motion judging means 24 
judges the start of the double stance state and the finish of the single stance state. 

[0070] Further, in the processing of the leg body motion judging means 24, of estimated values of floor reaction 

20 forces Ff and Fr (see Fig. 1(b)) acting on both leg bodies 2, 2, respectively, which are derived in the double stance 
state by the floor reaction force estimating means 35 as will be described later, the estimated value of the floor reaction 
force Fr= T (Frx, Frz) (specifically, an absolute value =V(Frx 2 +Frz 2 ) of the floor reaction force Fr derived at the previous 
cycle time of the arithmetic processing unit 1 6) acting on the leg body 2 on a rear side relative to the advancing direction 
of the human being is compared with a prescribed threshold value (positive value approximately "0") determined in 

25 advance. Then, when the absolute value of the estimated value of the floor reaction force Fr is lowered to the threshold 
value or less, it is judged that the double stance state finishes and the single stance state starts. In the present em- 
bodiment, the initial state of the motion state of the leg bodies 2, 2 is the double stance state, and the leg body motion 
judging means 24 judges that the motion state of the leg bodies 2, 2 is the double stance state until an estimated value 
of a floor reaction force acting on either one of the leg bodies 2, 2 is lowered to the foregoing threshold value or less. 

30 [0071] in parallel with the processing of the leg body motion judging means 24 as described above, the arithmetic 
processing unit 1 6 executes the processing by the foregoing chest inclination angle measuring means 25 and waist 
inclination angle measuring means 26. In this case, in the processing of the chest inclination angle measuring means 
25, an inclination angle 9a of the chest 4 in the absolute coordinate system Cf is derived in sequence per the foregoing 
cycle time according to a known technique using the processing of a so-called Kaiman filter, based on detection data 

35 about the acceleration of the chest 4 in the forward/backward direction and the angular velocity of the chest 4 which 
are inputted from the chest forward/backward accelerometer 15 and the chest gyro sensor 14, respectively. Likewise, 
in the processing of the waist inclination angle measuring means 25, an inclination angle Gb of the waist 3 in the absolute 
coordinate system Cf is derived in sequence using the processing of the Kaiman filter, based on detection data about 
the acceleration of the waist 3 in the forward/backward direction and the angular velocity of the waist 3 which are 

^0 inputted from the waist forward/backward accelerometer 20 and the waist gyro sensor 19, respectively. Here, in the 
present embodiment, the inclination angles Ga and Gb respectively of the chest 4 and the waist 3 in the absolute 
coordinate system Cf are inclination angles relative to, for example, the vertical direction (direction of gravity). 
[0072] It is also possible to derive inclination angles of the chest 4 and the waist 3 by, for example, integrating detection 
data about angular velocities by the gyro sensors 14 and 1 9. However, like in the present embodiment, the inclination 

45 angles Ga and Gb of the chest 4 and the waist 3 can be accurately measured using the processing of the Kaiman filter. 
[0073] Then, the arithmetic processing unit 1 6 executes the processing of the foregoing leg body posture calculating 
means 28 and the processing of the foregoing reference acceleration measuring means 27. 

[0074] tn the processing by the foregoing leg body posture calculating means 28, inclination angles Gc and 6d (in- 
clination angles relative to the vertical direction; see Fig. 2) of the thigh portion 9 and the crus portion 11 of each leg 
50 body 2 in the absolute coordinate system Cf are derived per the foregoing cycle time in the following manner. Specif- 
ically, the inclination angle Gc of the thigh portion 9 of each leg body 2 is calculated from a current-time value of detection 
data about the bending angle A9c of the hip joint 8 by the foregoing hip joint angle sensor 22 attached to the leg body 
2 and a current-time value of the inclination angle 0b of the waist 3 derived by the foregoing waist inclination angle 
measuring means 25, using the following equation (6). 



55 



Gc - Gb + AGc (6) 
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[0075] Here, the inclination angle 0b of the waist 3 takes a negative value when the waist 3 is inclined relative to the 
verticai direction such that an upper end portion of the waist 3 is projected forward of the human being 1 as compared 
with a lower end portion thereof, and the bending angle ABc of the hip joint 8 takes a positive value when the thigh 
portion 9 is inclined relative to the central axis of the waist 3 such that a lower end portion of the thigh portion 9 is 
5 projected forward of the human being 1 . 

[0076] Further the inclination angle Bd of Lhe crus portion 11 of each leg body 2 is calculated from a current-time 
value of the inclination angle Be of the thigh portion 9 derived in the foregoing manner and a current-lime value of 
detection data about the bending angle ABd of the knee joint 1 0 by the foregoing knee joint angle sensor 23 attached 
to the leg body 2. using the following equation (7). 

10 

Bd = Be - ABd (7) 

[0077] Here, the bending angle of the knee joint 1 0 takes a positive value when the crus portion 1 1 is inclined relative 

15 to the central axis of the thigh portion 9 toward the back side of the thigh portion 9. 

[0078] Further, in the processing of the foregoing reference acceleration measuring means 27, an acceleration a 0 = T 
(aOx, aOz) of the origin O of the foregoing body coordinate system Cp in the absolute coordinate system Cf is derived 
in the following manner. Specifically, given that a current-time value of detection data about the acceleration of the 
waist 3 in the forward/backward direction by the foregoing waist forward/backward accelerometer 20 is ap, and a 

20 current-time value of detection data about the acceleration of the waist 3 in the upward/downward direction by the 
foregoing waist upward/downward accelerometer 21 is aq, the acceleration a 0 = T (a 0 x, a 0 z) in the absolute coordinate 
system Cf is derived from those detection data ap and aq, and a current-time value of the inclination angle Bb of the 
waist 3 derived by the foregoing waist inclination angle measuring means 25, using the following equation (8). 



25 



a 0 = T (a 0 x, a 0 z) 



= T (ap • cosGb - aq • sinBb, ap • sinBb + aq • cosBb - g) (8) 

30 [0079] Then, the arithmetic processing unit 16 executes the processing of the foregoing respective-portion center of 
gravity position calculating means 29, and derives positions (positions relative to the origin of the body coordinate 
system Cp) of the centers of gravity of respective rigid-body corresponding portions of the human being 1 in the fore- 
going body coordinate system Cp, using a rigid body link model described hereinbelow. 

[0080] As shown in Fig. 4, a rigid body link model R used in the present embodiment is a model that expresses the 
35 human being 1 by linking together rigid bodies R1 , R1 corresponding to the thigh portions 9 of the respective leg bodies 
2. 2, rigid bodies R2, R2 corresponding to the crus portions 1 1 , 1 1 , a rigid body R3 corresponding to the waist 3, and 
a rigid body R4 corresponding to a portion 38 (hereinafter referred to as the upper body portion 38) in a combination 
of the foregoing chest 4, arm bodies 7, 7, and head 6. In this case, a linkage portion between each rigid body R1 and 
the rigid body R3, and a linkage portion between each rigid body R1 and the rigid body R2 correspond to the hip joint 
40 8 and the knee joint 10, respectively. Further, a linkage portion between the rigid body R3 and the rigid body R4 is a 
tilting fulcrum portion 39 for the chest 4 relative to the waist 3. 

[0081] In the present embodiment, positions of the centers of gravity G1 , G2, G3, and G4 of the rigid-body corre- 
sponding portions (the thigh portion 9 and the crus portion 11 of each leg body 2, the waist 3, and the upper body 
portion 38) of the human being 1 corresponding to the respective rigid bodies R1 to R4 of such a rigid body link model 
45 R, in the respective rigid-body corresponding portions, are derived in advance and stored in a memory, not shown, of 
the arithmetic processing unit 16. 

[0082] Here, the positions of the centers of gravity G1 , G2, G3, and G4 of the respective rigid-body corresponding 
portions stored/retained in the arithmetic processing unit 16 are positions each in a coordinate system that is fixed 
relative to each rigid-body corresponding portion . In this case, as data representing each of the positions of the centers 

50 of gravity G1 , G2, G3, and G4 of the respective rigid-body corresponding portions, for example, a distance in a direction 
of the central axis of the rigid-body corresponding portion from the central point of a joint at one end portion of such a 
rigid-body corresponding portion is used. Specifically, for example, as shown in Fig. 4, the position of the center of 
gravity G1 of each thigh portion 9 is given as a position with a distance t1 in the direction of the central axis of the thigh 
portion 9 from the center of the hip joint 8 of the thigh portion 9, the position of the center of gravity G2 of each crus 

55 portion 11 is given as a position with a distance t2 in the direction of the central axis of the crus portion 11 from the 
center of the knee joint 10 of the crus portion 11 , and values of those distances t1 and t2 are derived in advance and 
stored/retained in the arithmetic processing unit 1 6. This also applies to the positions of the centers of gravity G3, G4 
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of the other rigid-body corresponding portions. 

[0083] Strictly speaking, the position of the center of gravity G4 of the upper body portion 38 is affected by movement 
of the arm bodies 7 7 included in the upper body portion 38. However, inasmuch as the arm bodies 7, 7 during walking 
generally make a symmetrical positional relationship relative to the central axis of the chest 4, the position of the center 
of gravity G4 of the upper body portion 38 does not change much, and is, for example, substantially equivalent to a 
position of the center of gravity G4 of the upper body portion 38 in the erect posture stale. 

[0084] Further, in the present embodiment, in addition to the data representing the positions of the centers of gravity 
G1 , G2, G3, and G4 of the respective rigid-body corresponding portions (the thigh portion 9 and the crus portion 1 1 of 
each leg body 2, the waist 3, and the upper body portion 38), data about weights of the respective rigid-body corre- 
sponding portions, and data about sizes of the respective rigid-body corresponding portions (e.g. data about lengths 
of the respective rigid-body corresponding portions) are derived in advance and stored/retained in the arithmetic 
processing unit 1 6. 

[0085] The weight of the crus portion 11 is a weight including the foot portion 13. Further, the data stored/retained 
beforehand in the arithmetic processing unit 16 as described above may be obtained by actual measurement or the 
like, or may be estimated from a height and weight of the human being 1 based on human average statistical data In 
general, the positions of the centers of gravity G1 , G2,G3, and G4, the weights, andthe sizes of the foregoing respective 
rigid-body corresponding portions have correlation with a height or weight of a human being and, based on the corre- 
lation, it is possible to estimate relatively accurately the positions of the centers of gravity G1 , G2, G3 ; and G4 ( the 
weights, and the sizes of the foregoing respective rigid-body corresponding portions from data about the height and 
weight of the human being. 

[0086] The foregoing respective-portion center of gravity position calculating means 29 derives positions of the cent- 
ers of gravity G1 , G2, G3, and G4 of the respective rigid-body corresponding portions in the body coordinate system 
Cp (xz coordinates in Fig, 4) having the origin O fixed to the waist 3, from the data stored/retained beforehand in the 
arithmetic processing unit 16 as described above, current-time values of the inclination angle Ga of the chest 4 (=the 
inclination angle of the upper body portion 38) and the inclination angle 6b of the waist 3 respectively derived by the 
foregoing chest inclination angle measuring means 25 and waist inclination angle measuring means 26, and current- 
time values of the inclination angles Gc and Gd respectively of the thigh portion 9 and the crus portion 11 of each leg 
body 2 derived by the foregoing leg body posture calculating means 28. 

[0087] In this case, since the inclination angles Ga to Gd of the respective rigid-body corresponding portions (the thigh 
portion 9 and the crus portion 1 1 of each leg body 2. the waist 3, and the upper body portion 38) are derived as described 
above, the positions and postures of the respective rigid-body corresponding portions in the body coordinate system 
Cp can be derived from data about the inclination angles Ga to Gd and data about the sizes of the respective rigid-body 
corresponding portions. Therefore, the positions of the centers of gravity G1 , G2, G3, and G4 of the respective rigid- 
body corresponding portions in the body coordinate system Cp are derived. 

[0088] Specifically, for example, referring to Fig. 4, with respect to the leg body 2 located on the left side in Fig. 4, 
since an inclination angle of the thigh portion 9 in the body coordinate system Cp (an inclination angle relative to the 
z-axis direction) is 9c (in this case, Gc<0 in Fig. 4), coordinates of the position of the center of gravity G1 of the thigh 
portion 9 in the body coordinate system Cp become (t1 - sin Gc, -t1 • cos Gc). On the other hand, since an inclination 
angle of the crus portion 11 in the body coordinate system Cp is Gd (Gd<0 in Fig. 4), given that a length of the thigh 
portion 9 is Lc. coordinates of the position of the center of gravity G2 of the crus portion 11 in the body coordinate 
system Cp become (Lc - sin Gc +- t2 • sin Gd, -Lc • cos Gc - 12 • cos Gd). The centers of gravity of the thigh portion 9 and 
the crus portion 1 1 of the other leg body 2, the waist 3, and the upper body portion 38 can also be derived as described 
above. 

[0089] After deriving the positions of the centers of gravity G1 , G2, G3, and G4 of the respective rigid-body corre- 
sponding portions in the body coordinate system Cp by the respective-portion center of gravity position calculating 
means 29 as described above, the arithmetic processing unit 1 6 executes the processing of the foregoing body center 
of gravity position calculating means 30, and derives a position (xg, zg) of the body center of gravity GO of the human 
being 1 in the body coordinate system Cp using data about the positions of the centers of gravity G1 , G2, G3, and G4 
of the respective rigid-body corresponding portions and data about the weights of the respective rigid-body correspond- 
ing portions. 

[0090] Here, given that the position of the center of gravity G3 of the waist 3 in the body coordinate system Cp and 
the weight thereof are (x3, z3) and m3, respectively, the position of the center of gravity G4 of the upper body portion 
38 and the weight thereof are (x4, z4) and m4, respectively, the position of the center of gravity G1 of the thigh portion 
9 of the leg body 2 on the left side of the human being 1 facing forward and the weight thereof are (x1 L, z1L) and ml L, 
respectively, the position of the center of gravity G2 of the crus portion 1 1 of the same leg body 2 and the weight thereof 
are (x2L, z2L) and m2L, respectively, the position of the center of gravity G1 of the thigh portion 9 of the leg body 2 on 
the right side and the weight thereof are (x1 R, z1 R) and ml R, respectively, the position of the center of gravity G2 of 
the crus portion 1 1 of the same leg body 2 and the weight thereof are (x2R, z2R) and m2R, respectively and the weight 
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of the human being 1 is M (=m1 L+m2L+m1 R+m2R+m3+m4), the position (xg } zg) of the body center of gravity GO of 
the human being 1 in the body coordinate system Cp is derived by the following equations (9). 

5 xg = (ml L • x1 L + ml R . x1 R + m2L ■ x2L + m2R . x2R 

+• m3 • x3 + m4 • x4) / M 
zg = (ml L • z1 L + ml R • z1 R + m2L ■ z2L + m2R ■ z2R 

+ m3-z3 + m4-z4)/M (9) 

[0091 ] After executing the processing of the body center of gravity position calculating means 30 as described above, 
the arithmetic processing unit 16 further executes the processing of the foregoing body center of gravity acceleration 

15 calculating means 32 and the processing of the foregoing ankle position calculating means 31 . 

[0092] In this case, in the processing of the body center of gravity acceleration calculating means 32, at the outset, 
using time series data about the positions (xg, zg) of the body center of gravity GO in the body coordinate system Cp 
derived by the body center of gravity position calculating means 30 per the foregoing cycle time, a second-order dif- 
ferentiated value of the position (xg, zg) of the body center of gravity GO in the body coordinate system Cp, that is, an 

20 acceleration T (d 2 xg/dt 2 , d 2 zg/dt 2 ) of the body center of gravity GO relative to the origin O of the body coordinate system 
Cp, is derived. Then, by deriving the vector sum of the acceleration T (d 2 xg/dt 2 , d 2 zg/dt 2 ) and the acceleration a 0 = T 
(a 0 x, a 0 z), derived by the foregoing reference acceleration measuring means 27, of the origin O of the body coordinate 
system Cp in the absolute coordinate system Cf , an acceleration a= T (ax, az) of the body center of gravity GO in the 
absolute coordinate system Cf is derived. 

25 [0093] Further, in the processing of the foregoing ankle position calculating means 31 , at the outset, from current- 
time values of data about the inclination angles Be and 6d respectively of the thigh portion 9 and the cms portion 1 1 of 
each leg body 2 derived by the foregoing leg body posture calculating means 28, a current-time value of data about 
the inclination angle 6b of the waist 3 derived by the foregoing waist inclination angle measuring means 26, and data 
about the sizes (lengths) of the thigh portion 9 and the crus portion 11 , a position of the ankle portion 12 of each leg 

30 body 2 in the foregoing body coordinate system Cp is derived by the processing like the processing of the foregoing 
respective-portion center of gravity position calculating means 29. Specifically, referring to Fig. 4, with respect to the 
leg body 2 located on the left side in Fig. 4, given that a length of the crus portion 11 (a length from the center of the 
knee joint 10 to the ankle portion 12) is Ld, coordinates (x12, z12) of the position of the ankle portion 12 in the body 
coordinate system Cp become (Lc • sin 6c + Ld • sin 0d, - Lc • cos 6c - Ld - cos 6d) (wherein 6c<0 and 6d<0 in Fig. 4). 

35 This also applies to the other leg body 2. 

[0094] Then, from current-time values of data about the position (x12, z12) of the ankle portion 12 in the body coor- 
dinate system Cp and the position (xg, zg) of the body center of gravity GO in the body coordinate system Cp derived 
by the foregoing body center of gravity position calculating means 30, a position vector T (x12-xg, z12-zg) of the ankle 
portion 12 of each leg body 2 relative to the body center of gravity GO, that is, AXf, AZf, AXr, and AZr in the foregoing 

40 equations (5), are derived. 

[0095] Then, the arithmetic processing means 1 6 executes the processing of the foregoing floor reaction force es- 
timating means 36 in the following manner. Specifically, in this processing, when the motion state of the leg bodies 2, 
2 determined by the foregoing leg body motion judging means 24 at the current cycle time is the single stance state, 
an estimated value of a floor reaction force F= T (Fx, Fz) acting on the leg body 2 touching the ground is derived from 

45 values of the weight M and the acceleration of gravity g of the human being 1 (these are stored beforehand in the 
arithmetic processing unit 1 6) and a current-time value of the acceleration a= T (ax, az) of the body center of gravity GO 
in the absolute coordinate system Cf derived by the foregoing body center of gravity acceleration calculating means 
32, using the foregoing equation (2). In this case, a floor reaction force acting on the ground non-touching side leg 
body 2 (the free-leg side leg body 2) is T (0, 0). 

50 [0096] On the other hand, when the motion state of the leg bodies 2. 2 determined by the leg body motion judging 
means 24 at the current cycle time is the double stance state, estimated values of floor reaction forces Ff= T (Ffx, Ffz) 
and Fr= T (Frx, Frz) on the respective leg bodies 2, 2 are derived from the weight M and the acceleration of gravity g of 
the human being 1 , a current-time value of the acceleration a= T (ax, az) of the body center of gravity GO in the absolute 
coordinate system Cf derived by the foregoing body center of gravity acceleration calculating means 32, and current- 

55 time value data about the position of the ankle portion 12 of each leg body 2 relative to the body center of gravity GO 
(current-time values of data about AXf, AZf, AXr, and AZr in the equations (5)) derived by the foregoing ankle position 
calculating means 31 , using the foregoing equations (5). 

[0097] On the other hand, in parallel to the processing of the body center of gravity position calculating means 30, 
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the body center of gravity acceleration calculating means 32, the ankle position calculating means 31 , and the floor 
reaction force estimating means 36 as described above, the arithmetic processing unit 16 executes the processing of 
the foregoing leg body respective-portion acceleration calculating means 33, leg body respective-portion angular ve- 
locity calculating means 34, and floor reaction force acting point estimating means 35. 

5 [0098] In this case, in the processing of the foregoing leg body respective-portion acceleration calculating means 
33, al the outset, like in the foregoing body center of gravity acceleration calculating means 32, using time series data 
about the positions of the centers of gravity G1 and G2 respectively of the thigh portion 9 and the crus portion 1 1 being 
the rigid-body corresponding portions of each leg body 2 in the body coordinate system Cp, which are derived by the 
foregoing respective-portion center of gravity position calculating means 29 per the foregoing cycle time, second-order 

10 differentiated values of the positions of the centers of gravity G1 and G2 respectively of the thigh portion 9 and the 
crus portion 11 in the body coordinate system Cp, that is, accelerations (accelerations relative to the origin O of the 
body coordinate system Cp) of the respective centers of gravity G1 and G2 of the thigh portion 9 and the crus portion 
1 1 in the body coordinate system Cp, are derived. Then by deriving the vector sum of these respective accelerations 
and the acceleration a 0 -^ T (a 0 x, a 0 z) of the waist 3 in the absolute coordinate system Cf by the foregoing reference 

15 acceleration measuring means 27, accelerations (more specifically, coordinate components of the accelerations in the 
absolute coordinate system Cf) respectively of the thigh portion 9 and the crus portion 11 in the absolute coordinate 
system Cf are derived. 

[0099] Further, in the processing of the foregoing leg body respective-portion angular velocity calculating means 34, 
using time series data about the inclination angles 9c and Bd respectively of the thigh portion 9 and the crus portion 

20 11 of each leg body 2 derived by the foregoing leg body posture calculating means 28 per the foregoing cycle time, 
second-order differentiated values of the respective inclination angles Be and Bd of the thigh portion 9 and the cms 
portion 1 1 , that is, angular velocities respectively of the thigh portion 9 and the crus portion 11 , are derived. 
[0100] Further, in the processing of the floor reaction force acting point estimating means 35, with respect to the leg 
body 2 touching the ground, a vector (position vector of a floor reaction force acting point relative to the ankle portion 

25 12; hereinafter referred to as the floor reaction force acting point vector) from the ankle portion 12 of the leg body 2 to 
a floor reaction force acting point of the foot portion 13 of the leg body 2 (a point upon which the total floor reaction 
force acting on a ground touching portion of the foot portion 13 is considered to be concentrated) is derived as data 
representing a position of the floor reaction force acting point, from, for example, a current-time value of the inclination 
angle 6c of the thigh portion 9 derived by the foregoing leg body posture calculating means 28, based on predetermined 

30 correlations as shown in Figs. 5 and 6. 

[0101] Specifically, according to knowledge of the inventors of the present application, the inclination angle 6c of the 
thigh portion 9 or the bending angle ABd of the knee joint 10 of the leg body 2 touching the ground has relatively 
significant correlation with the floor reaction force acting point. The foregoing floor reaction force acting point vector, 
specifically, components of the floor reaction force acting point vector in the advancing direction (X-axis direction) of 

35 the human being 1 , and components of the floor reaction force acting point vector in the vertical direction (Z-axis vector), 
respectively change with respect to, for example, the inclination angle ec of the thigh portion 9 as shown in FIGS. 5 
and 6. Here, the negative inclination angle 6c of the thigh portion 9 is an angle when the thigh portion 9 is inclined 
relative to the central axis of the waist 3 so that the leg body 2 extends to the back side of the human being 1 (e.g. the 
leg body 2 on the right side of the human being 1 facing forward in Fig. 2), while the positive inclination angle Gc is an 

40 angle when the thigh portion 9 is inclined relative to the central axis of the waist 3 so that the leg body 2 exists on the 
front side of the human being 1 (e.g. the leg body 2 on the left side of the human being 1 facing forward in Fig. 2). 
[0102] Therefore, in the present embodiment, approximate expressions using the inclination angle 0c of the thigh 
portion 9 as a parameter, which represent the correlations of Figs. 5 and 6, are prepared and stored/retained in advance 
in the arithmetic processing unit 1 6. Then, in the processing of the foregoing floor reaction force acting point estimating 

45 means 35, a current-time value of the inclination angle Be of the thigh portion 9 derived by the foregoing leg body 
posture calculating means 28 is put in the foregoing approximate expressions to thereby derive the foregoing floor 
reaction force acting point vector (specifically, components of the floor reaction force acting point vector in the X-axis 
direction and Z-axis direction). 

[0103] Here, in the correlations as shown in Figs. 5 and 6 wherein the inclination angle ec of the thigh portion 9 has 
50 a minimal value, even when the inclination angle Be of the thigh portion 9 is the same, a value of the floor reaction 
force acting point vector differs between a decreasing stage of the inclination angle Be and an increasing stage thereof. 
Therefore, in the present embodiment, upon preparing the foregoing approximate expression, a transition of the fore- 
going correlation from a landing of the heel of the foot portion 1 3 on the floor to a takeoff of the tiptoe from the floor is 
classified into a first phase (a phase of a1 in Fig. 5; a phase of b1 in Fig. 6) wherein the inclination angle Be of the thigh 
55 portion 9 is positive, a second phase (a phase of a2 in Fig, 5; a phase of b2 in Fig. 6) wherein the inclination angle Be 
of the thigh portion 9 is negative, and a changing speed of the inclination angle Be of the thigh portion 9, i.e. an inclination 
angular velocity of the thigh portion 9, is negative, and a third phase (a phase of a3 in Fig. 5; a phase of b3 in Fig. 6) 
wherein the inclination angle 9c of the thigh portion 9 is negative, and an Inclination angular velocity of the thigh portion 
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9 is positive, and the phases are approximated by the same function or different functions with respect to X-axis direction 
components and Z-axis direction components of the floor reaction force acting point vector, respectively. An approxi- 
mate expression of a phase in a combination of the first and second phases a1 and a2 in the correlation of Fig. 5 is, 
given that the X-axis direction components of the floor reaction force acting point vector are px, expressed by ; for 
example, a sixth-order polynomial function in the form of 

px = x 1 • 9c 6 + x 2 • 9c 5 + x 3 • Be 4 

3 2 

+ x 4 • 6c 4- x 5 • 6c + x 6 • Be 4- x 7 

(x 1 to x 7 are constant values). 

Further, an approximate expression of the third phase a3 in the correlation of Fig. 5 is expressed by, for example, a 
fourth-order polynomial function in the form of 

p x = x 8 • Be 4 + x 9 - 6c° + x 10 - Be 2 + x 1 1 • Be 4- x 12 
(x 8 to x 12 are constant values). 

20 [0104] Further, an approximate expression of a phase in a combination of the first and second phases b1 and b2 in 
the correlation of Fig. 6 is, given that the Z-axis direction components of the floor reaction force acting point vector are 
pz, expressed by, for example, a sixth-order polynomial function in the form of 



10 



15 



25 



pz = z 1 • Be 6 4- z 2 • Be 5 4- z 3 • Be 4 



4- Z A • BC 3 4- Z R • BC 2 4- Z R • BC 4- Z 7 



(z 1 to z 7 are constant values). 

30 Further, an approximate expression of the third phase b3 in the correlation of Fig. 6 is expressed by, for example, a 
third-order polynomial function in the form of 

3 2 

pz = z 8 ■ Be + z 9 ■ 0c ■ + z 10 • 6c 4- z^ 

35 

(x 8 to x^ are constant values). 
[0105] Then, upon deriving the floor reaction force acting point vector it is judged whether the inclination angle Be 
of the thigh portion 9 is positive or negative, and further judged whether the inclination angular velocity of the thigh 
portion 9 calculated by first-order differentiation of time series data about the inclination angles 6c of the thigh portion 

40 9 is positive or negative. Further, from judged positive or negative of the inclination angle 6c and judged positive or 
negative of the inclination angular velocity, it is judged which of the phases is manifested currently and, by putting a 
current-time value of the inclination angle 9c of the thigh portion 9 into the approximate expression of the judged phase, 
the floor reaction force acting point vector is calculated. Thereby, the value of the floor reaction force acting point vector 
in the decreasing stage of the inclination angle 6c of the thigh portion 9, and the value of the floor reaction force acting 

45 point vector in the increasing stage thereof can be discriminatingly calculated. 

[0106] In the present embodiment, the correlation between the inclination angle 6c of the thigh portion 9 of the leg 
body and the floor reaction force acting point vector is approximated by the polynomials to thereby derive the floor 
reaction force acting point vector. However, it is also possible to store/retain the correlations shown in Figs. 5 and 6 in 
the form of a data table and, using the data table, derive the floor reaction force acting point vector from the inclination 

so angle 6c of the thigh portion 9. 

[0107] Further, the position of the floor reaction force acting point also has correlation with the bending angle of the 
knee joint 1 0 of the leg body 2 touching the ground. Therefore, the position of the floor reaction force acting point may 
be estimated from the bending angle A6d of the knee joint 1 0 measured by the knee joint angle sensor 23, instead of 
the inclination angle Be of the thigh portion 9. Alternatively, using both the inclination angle Be of the thigh portion 9 

55 and the bending angle ABd of the knee joint 1 0, the position of the floor reaction force acting point may be estimated 
by the use of a map or the like. 

[0108] Further, when the human being 1 sits in a chair or stands up from the state of sitting in the chair, a correlation 
shown in Fig. 7 (upon sitting in the chair) or Fig. 8 (upon standing up from the chair) is established between the position 
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of the floor reaction force acting point and the bending angle AOd of the knee joint 10 and, when going upstairs or 
downstairs. Fig. 9 (upon going upstairs) or Fig. 1 0 (upon going downstairs) is established between the position of the 
floor reaction force acting point and the inclination angle Be of the thigh portion 9. Therefore, upon sitting in the chair 
or standing up ; the position of the floor reaction force acting point can be estimated from the bending angle AB6 of the 
5 knee joint 1 0 based on the correlation of Fig. 7 or Fig. 8, while, upon going upstairs or downstairs, the position of the 
floor reaction force acting point can be estimaled from the inclination angle 6c of the thigh portion 9 based on the 
correlation of Fig. 9 or Fig. 10. 

[0109] After estimating the position of the floor reaction force acting point as described above, then, the arithmetic 
processing unit 1 6 executes the processing of the foregoing joint moment estimating means 37, and derives moments 

10 acting on the knee joint 10 and the hip joint 8 of each leg body 2. This processing is carried out based on a so-called 
inverse dynamics model using current-time values of the data respectively derived by the foregoing floor reaction force 
estimating means 36, leg body respective-portion acceleration calculating means 33, leg body respective-portion an- 
gular velocity calculating means 34, floor reaction force acting point estimating means 35, and leg body posture cal- 
culating means 28. This inverse dynamics model uses an equation of motion for translation motion of each rigid-body 

15 corresponding portion of the human being 1 and an equation of motion for rotating motion thereof so as to derive 
moments acting on joints in order from the joint closer to the floor reaction force acting point. In the present embodiment, 
the moments acting on the knee joint 10 and the hip joint 8 of each leg body 2 are derived in order. 
[0110] More specifically, referring to Fig. 11 , at the outset, with respect to the cms portion 11 of each leg body 2, a 
force (joint reaction force) acting on the ankle portion 12 at a tip portion of the cms portion 11 f a force (joint reaction 

20 force) acting on the portion of the knee joint 1 0 of the crus portion 1 1 , and a translation acceleration of the center of 
gravity G2 of the crus portion 11 are respectively given as T (F-,x, FiZ), T (F 2 x, F 2 z) : and T (a 2 x, a 2 z) according to the 
component notation in the absolute coordinate system Cf, and the weight of the crus portion 1 1 is given as m 2 . In this 
case, an equation of motion for the translation motion of the center of gravity G2 of the crus portion 11 becomes the 
following equation (10). 

25 

T (m 2 • a 2 x, m 2 • a 2 z) ~ T (F 1 x-F 2 x, F 1 z-F 2 z - m 2 • g) 

therefore, 

30 

T (F 2 x, F 2 z) = T (F 1 x - m 2 • a 2 x, F-,z - m 2 • a 2 z - m 2 • g) (10) 

[01 11 ] Here, the acceleration T (a 2 x, a 2 z) of the center of gravity G2 of the crus portion 1 1 is derived by the foregoing 
3 5 leg body respective-portion acceleration calculating means 33. Further, the joint reaction force T (F-,x, F-|Z) acting on 
the ankle portion 12 at the tip portion of the crus portion 11 is approximately equal to the estimated value of the floor 
reaction force derived by the foregoing floor reaction force estimating means 36 for the leg body 2 having the crus 
portion 11 . More specifically, in the single stance state, when the leg body 2 is on the ground, the joint reaction force T 
(F^, F^) is the floor reaction force T (Fx, Fz) derived by the foregoing equation (2) and, when the leg body 2 is the 
^0 f ree _| e g side | e g body, T (F 1 x, F 1 z)= T (0, 0). On the other hand, in the double stance state, when the leg body 2 is the 
leg body on the rear side of the human being 1 facing forward in the advancing direction, the joint reaction force T (F-,x, 
F 1 z) is the floor reaction force T (Frx, Frz) of the foregoing equations (5) and, when the leg body 2 is the front-side leg 
body, it is the floor reaction force T (Ffx, Ffz) of the foregoing equations (5). 

[0112] Therefore, the joint reaction force T (F 2 x, F 2 z) acting on the knee joint 10 of each leg body 2 is derived from 
45 data about the acceleration T (a 2 x, a 2 z) of the center of gravity G2 of the crus portion 11 derived by the leg body re- 
spective-portion acceleration calculating means 33, data about the floor reaction force (= T (F-,x, F 1 z)) derived by the 
floor reaction force estimating means 36, data about the weight m 2 of the crus portion 11 derived in advance, and a 
value of the acceleration of gravity g, using the foregoing equation (10). 

[0113] Further, referring to Fig. 11, it is given that a moment acting on the ankle portion 12 at the tip portion of the 
50 cms portion 11 is M.,, a moment acting on the portion of the knee joint 10 of the crus portion 11 is M 2: a moment of 
inertia about the center of gravity G2 of the crus portion 11 is l G2 , and an angular velocity about the center of gravity 
G2 of the crus portion 11 is a 2 Further, given that, in correspondence with the foregoing figure 4, a distance between 
the center of gravity G2 of the crus portion 11 and the center of the knee joint 10 is t 2 , and a distance between the 
center of gravity G2 of the crus portion 1 1 and the ankle portion 1 2 is t 2 * (=Ld-t 2 ), an equation of motion for the rotating 
55 motion about the center of gravity G2 of the crus portion 11 becomes the following equation (11). 

| Q2 . a 2 = M 1 - M 2 + F.,x • t 2 ' • cos ed - F.,z • t 2 ' * sin Bd 
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4- F 2 x ■ t 2 • cos Bd - F 2 7 • t 2 - sin Bd 



therefore 



M 2 = M 1 - l G2 • a 2 + F.,x ■ t 2 ' • cos Bd - F 1 z - 1 2 ' • sin Bd 



+ F 2 x • t 2 • cos Bd - F 2 z • t 2 - sin Bd 



(11) 



[0114] Here, M 1 in the equation (11) is a moment obtained as the outer product (vector product) of a floor reaction 
force acting point vector derived by the foregoing floor reaction force acting point estimating means 35 with respect to 
the leg body 2 having the crus portion 11 relating to the equation (11), and a floor reaction force vector derived by the 
foregoing floor reaction force estimating means 36 with respect to that leg body 2. ot 2 is an angular velocity of the crus 
portion 11 derived by the foregoing leg body respective-portion angular velocity calculating means 34. Bd is an incli- 
nation angle of the crus portion 1 1 derived by the foregoing leg body posture calculating means 28. T (F-,x, F^) is, as 
described before, the estimated value of the floor reaction force derived by the floor reaction force estimating means 
36. Further, T (F 2 x, F 2 z) is derived by the foregoing equation (10). The moment of inertia l G2 is derived in advance and 
stored in the arithmetic processing unit 1 6 along with data about the weight m 2 and size of the crus portion 1 1 , and so 
forth. 

[01 1 5] Therefore, the moment M 2 acting on the knee joint 1 0 is derived from data about the estimated value of the 
floor reaction force by the floor reaction force estimating means 36, data about the floor reaction force acting point 
vector by the floor reaction force acting point estimating means 35, data about the angular velocity ot 2 of the crus portion 
1 1 by the leg body respective-portion angular velocity calculating means 34, data about the inclination angle Bd of the 
crus portion 11 by the leg body posture calculating means 28, data about the joint reaction force T (F 2 x, F 2 z) derived 
by the foregoing equation (10), and data about the moment of inertia l G2; the size (Ld), and the position (t 2 ) of the 
center of gravity G2 of the crus portion 11 derived in advance, using the foregoing equation (11). 
[0116] After deriving the moment M 2 acting on the portion of the knee joint 10 of the crus portion 11 as described 
above, the joint moment estimating means 37 derives, according to the processing like the processing of calculation 
thereof, a moment acting on the portion of the hip joint 8 of the thigh portion 9. The basic idea of this processing is the 
same as the technique of deriving the moment M 2 of the knee joint 1 0, and therefore, detailed showing and explanation 
are omitted. An outline thereof is as follows. 

[0117] Specifically, at the outset, a joint reaction force T (F 3 x, F 3 z) acting on the portion of the hip joint 8 of the thigh 
portion 9 is derived using the following equation (12) (an equation in the same form as the foregoing equation (10)) 
based on an equation of motion for the translation motion of the center of gravity G1 (see Fig. 4) of the thigh portion 9. 



[0118] Here, T (F 2 x, F 2 z) is the joint reaction force of the knee joint 10 derived before by the foregoing equation (10). T 
(a^, a A z) is an acceleration (translation acceleration) of the center of gravity G1 of the thigh portion 9 in the absolute 
coordinate system Cf, which is derived by the foregoing leg body respective-portion acceleration calculating means 
33. Further, m-, is the weight of the thigh portion 9 derived in advance, and g is the acceleration of gravity. 
[0119] Then, a moment M 3 acting on the portion of the hip joint 8 of the thigh portion 9 is derived using the following 
equation (13) (an equation in the same form as the foregoing equation (11)) based on an equation of motion for the 
rotating motion about the center of gravity G1 of the thigh portion 9. 



[0120] Here, M 2 is the moment on the knee joint 1Q derived by the foregoing equation (11), T (F 2 x, F 2 z) is the joint 
reaction force on the knee joint 10 derived by the foregoing equation (10), T (F 3 x, F 3 z) is the joint reaction force on the 
hip joint 8 derived by the foregoing equation (12), l G1 is a moment of inertia about the center of gravity G1 of the thigh 
portion 9 derived in advance, ct-i is an angular velocity of the thigh portion 9 derived by the foregoing leg body respective- 
portion angular velocity calculating means 34, and Be is an inclination angle of the thigh portion 9 derived by the 



T (F 3 x. F 3 z) = T (F 2 x - m 1 ■ a.,x. F 2 z - m 1 ■ a A z - m 1 • g) 



(12) 



M 3 = M 2 - l G1 • a 1 + F 2 x • t-,' • cos 6c - F 2 z ■ t-,' • sin Be 



+ F 3 x • t 1 • cos Be - F 3 z • t 1 ■ sin Be 



(13) 
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foregoing leg body posture calculating means 28. Further t-, is a distance from the center of the hip joint 8 to the center 
of gravity G1 of the thigh portion 9 (see Fig. 4), and V is a distance from the center of the knee joint 10 to the center 
of gravity G1 of the thigh portion 9 (Lc-t 1 in Fig. 4), which are determined based on the position of the center of gravity 
G1 and the size (length) of the thigh portion 9 derived in advance. 
5 [0121] The processing described above is successively executed per cycle time of the foregoing arithmetic process- 
ing unit 1 6, so that the floor reaction force acting on each leg body Z, and the moments acting on the knee joint 10 and 
the hip joint 8 of each leg body 2 are sequentially estimated in real time. 

[0122] Although detailed explanation is omitted in this specification, the derived estimated values of the moments 
on the knee joints 10 and the hip joints 8 are used for a control of , for example, a device assisting walking of the human 
10 being 1 (a device including electric motors or the like that can give auxiliary torques to the knee joints 1 0 and the hip 
joints 8). 

[0123] The states of temporal variations of estimated values of ftoor reaction forces (specifically, absolute values of 
the estimated values of the floor reaction forces) derived by the foregoing processing of the arithmetic processing unit 
16 are illustrated in Figs. 12 to 14 by solid lines. Further the states of temporal variations of estimated values of 

15 moments on the knee joint 10 and the hip joint 8 derived by the processing of the arithmetic processing unit 16 are 
illustrated in Fig. 15 by solid lines. Here, Figs. 12 and 15 are illustration wherein the human being 1 walks on a fiatland 
at a substantially constant speed, Fig. 13 is illustration wherein the human being 1 walks upstairs, and Fig. 14 is 
illustration wherein the human being 1 stands up from the state of sitting in a chair. In Figs. 12 to 14, comparative 
examples (corresponding to true values of floor reaction forces) obtained by actually measuring floor reaction forces 

20 using a force meter or the like are also shown by imaginary lines. Further, in Fig. 15, comparative examples (corre- 
sponding to true values of moments on the knee joint 1 0 and the hip joint 8) obtained by actually measuring moments 
on the knee joint 1 0 and the hip joint 8 using a torque meter or the like are also shown by imaginary lines. 
[0124] As is clear by referring to Figs. 12 to 14, it is understood that, according to the present embodiment, the 
accurate estimated values of the floor reaction forces are obtained irrespective of the motion type or motion environment 

25 of the leg bodies 2. Further, in the present embodiment, as shown in Fig. 15, the moments on the knee joint 10 and 
the hipjoint8can also be estimated with relatively high accuracy using the estimated values of the floor reaction forces. 
[0125] As described above, according to this embodiment, it is possible to easily estimate in real time the floor reaction 
force acting on each leg body 2 and the moments acting on the hip joint 8 and the knee joint 10 of each leg body 2 
using relatively small and light sensors such as the angle sensors 22 and 23 respectively attached to the hip joints 8 

30 and the knee joints 10, and the gyro sensors 14 and 19 and the accelerometers 15. 20, and 21 attached to the body 
5, without attaching such sensors to the leg bodies 2 that become obstructive to the walking of the human being 1 , or 
apply loads to the motion of the leg bodies'2. Moreover, such estimation can be carried out with relatively high accuracy 
irrespective of the motion type or motion environment of the leg bodies 2, i.e. whether it is walking on a fiatland, walking 
on the stairs, or the like. 

35 [0126] In the embodiment described above, the description has been given about the case, as an example, wherein 
the present invention is applied to the human being 1. However, the present invention is also applicable to a bipedal 
walking robot as a bipedal walking moving body. Here, in case of the bipedal walking robot, it may have a structure 
where a waist and a chest are integral with each other In this case, it is possible to attach the gyro sensor and the 
forward/backward accelerometer to either one of the waist and the chest, and estimate floor reaction forces and mo- 
ments on joints of leg bodies in the same manner as in the present embodiment. Further, in case of the bipedal walking 
robot, it is also possible grasp bending angles of hip joints and knee joints based on control amounts of a controller 
for actuators of those joints. 

[0127] In the foregoing embodiment, the detection data of the waist upward/downward accelerometer 21 is uses as 
it is forjudging the motion state of the leg bodies 2. However, instead of such detection data, it is possible to use, for 
45 example, values of components of the acceleration a 0 of the waist 3 in the vertical direction (Z-axis direction) in the 
absolute coordinate system Cf , which are derived by the foregoing reference acceleration measuring means 27. 

INDUSTRIAL APPLICABILITY 

50 [0128] As described above, the present invention is useful in that, when performing walking assistance or walking 
control of a bipedal walking moving body such as a human being or a robot, floor reaction forces and joint moments 
acting on leg bodies of the moving body can be grasped. 



55 Claims 

1. A method of estimating a floor reaction force acting on each of leg bodies of the bipedal walking moving body, 
comprising: 
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a first step of judging whether a motion state of the leg bodies of said bipedal walking moving body is a single 
stance state in which only one of the leg bodies touches the ground : or a double stance state in which both 
leg bodies touch the ground; 

a second step of sequentially deriving positions of the center of gravity of said bipedal walking moving body, 
and sequentially deriving accelerations of said center of gravity in an absolute coordinate system fixed relative 
to Lhe ground by the use of time series data about the positions of said center of gravity; 
a third step of sequentially deriving positions of a specific portion relative to said center of gravity at least in 
said double stance state, said specific portion predetermined in the neighborhood of a lower end portion of 
each leg body; 

a step of deriving estimated values of said floor reaction force acting on the leg body touching the ground 
sequentially, in the single stance state of said bipedal walking moving body, based on an equation of motion 
for said center of gravity expressed by a weight and an acceleration of gravity of the bipedal walking moving 
body, the acceleration of said center of gravity, and said floor reaction force acting on the leg body touching 
the ground; and 

a step of deriving estimated values of said floor reaction forces respectively acting on both leg bodies sequen- 
tially, in the double stance state of said bipedal walking moving body, based on an equation of motion for said 
center of gravity expressed by a weight and an acceleration of gravity of the bipedal walking moving body, the 
acceleration of said center of gravity, and said floor reaction forces respectively acting on both leg bodies, and 
an expression of relation between the position of said specific portion of each leg body relative to said center 
of gravity and said floor reaction force acting on said leg body, said expression of relation being determined 
based on assumption that said floor reaction force acting on each leg body acts from said specific portion of 
said leg body toward said center of gravity. 

A method according to claim 1 , wherein said specific portion of each leg body is an ankle portion of said leg body. 

A method according to claim 1 , further comprising a step of measuring an acceleration, in an upward/downward 
direction, of a lower portion of a body supported on both leg bodies via a hip joint of each leg body, said lower 
portion of the body located near said hip joints, 

wherein, in said first step, the motion state of said bipedal walking moving body is judged such that when 
the acceleration of said lower portion of the body in the upward/downward direction increases to a predetermined 
threshold value or more, said double stance state starts while said single stance state finishes and, when the 
estimated value of said floor reaction force acting on the leg body which precedingly makes a takeoff is lowered 
to a predetermined threshold value or less in said double stance state, said double stance state finishes while said 
single stance state starts. 

A method according to claim 1 , further comprising a step of respectively measuring an inclination angle of a body 
supported on both leg bodies via a hip joint of each leg body, bending angles respectively of at least the hip joint 
and a knee joint of each leg body, and an acceleration of a predetermined reference point of said bipedal walking 
moving body in said absolute coordinate system, 

wherein, in said second step, based on the inclination angle of said body, the bending angles respectively 
of said hip joints and said knee joints, a rigid body link model formed by expressing said bipedal walking moving 
body as a linked body of a plurality of rigid bodies, prederived weights of respective rigid-body corresponding 
portions of the bipedal walking moving body corresponding to the respective rigid bodies of said rigid body link 
model, and positions of the prederived centers of gravity of the rigid-body corresponding portions in the respective 
rigid-body corresponding portions, positions of the center of gravity of said bipedal walking moving body relative 
to said reference point are sequentially derived, accelerations of said center of gravity relative to said reference 
point are sequentially derived based on time series data about the positions of said center of gravity, and the 
acceleration of said center of gravity in said absolute coordinate system is derived from the acceleration of said 
center of gravity relative to said reference point, and the acceleration of said reference point in said absolute 
coordinate system. 

A method according to claim 4, wherein said reference point is set to said body. 

A method according to claim 4, wherein said body has a waist coupled to both leg bodies via the hip joints, and a 
chest located on said waist so as to be tiltable relative to said waist, and the inclination angle of said body used 
for deriving the position of said center of gravity comprises inclination angles respectively of said waist and said 
chest. 
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A method according to claim 6 ; wherein said rigid body link model is a model expressing a crus portion located on 
a lower side of the knee joint of each leg body of said bipedal walking moving body ; a thigh portion between said 
knee joint and said hip joint, said waist, and an upper body portion located on an upper side of said waist and 
including said chest, as rigid bodies, respectively. 

A method of estimating a moment acting on at least one joint of each leg body of the bipedal walking moving body 
by the use of the estimated values of the floor reaction force on each leg body sequentially derived by the method 
according to claim 1 , comprising: 

a step of respectively measuring an inclination angle of a body supported on both leg bodies via a hip joint of 
each leg body, bending angles respectively of at least the hip joint and a knee joint of each leg body, and an 
acceleration of a predetermined reference point of said bipedal walking moving body in said absolute coordi- 
nate system; 

a step of sequentially deriving inclination angles of respective rigid-body corresponding portions of the bipedal 
walking moving body corresponding to the respective rigid bodies of said rigid body Sink model, based on the 
inclination angle of said body, the bending angles respectively of said hip joint and said knee joint of each leg 
body, and a rigid body link model formed by expressing said bipedal walking moving body as a linked body of 
a plurality of rigid bodies; 

a step of sequentially deriving positions of the centers of gravity of the respective rigid-body corresponding 
portions relative to said reference point, based on the inclination angles of said respective rigid-body corre- 
sponding portions, prederived weights of the respective rigid-body corresponding portions, and positions of 
the prederived centers of gravity of the rigid-body corresponding portions in the respective rigid-body corre- 
sponding portions, and sequentially deriving accelerations of the centers of gravity of the respective rigid-body 
corresponding portions relative to said reference point based on time series data about the positions of the 
centers of gravity of the respective rigid-body corresponding portions; 

a step of sequentially deriving accelerations of the centers of gravity of the respective rigid-body corresponding 
portions in said absolute coordinate system from the accelerations of the centers of gravity of the respective 
rigid-body corresponding portions relative to said reference point, and the accelerations of said reference point 
in said absolute coordinate system; 

a step of sequentially deriving angular velocities of the respective rigid-body corresponding portions based on 
time series data about the inclination angles of the respective rigid-body corresponding portions; 
a step of sequentially deriving estimated positions of a floor reaction force acting point of each leg body in 
said bipedal walking moving body, based on at least one of an inclination angle of a thigh portion of said leg 
body and the bending angle of the knee joint of said leg body as the rigid-body corresponding portions of said 
bipedal walking moving body; and 

a step of estimating a moment acting on at least one of the joints of each leg body of said bipedal walking 
moving body based on an inverse dynamics model using the estimated value of said floor reaction force, the 
estimated position of said floor reaction force acting point, the accelerations of the centers of gravity of the 
respective rigid-body corresponding portions and the angular velocities of said rigid-body corresponding por- 
tions in said absolute coordinate system, the inclination angles of the respective rigid-body corresponding 
portions, the prederived weights and sizes of the respective rigid-body corresponding portions, the positions 
of the prederived centers of gravity of the rigid-body corresponding portions in the respective rigid-body cor- 
responding portions, and prederived moments of inertia of the respective rigid-body corresponding portions. 

A method of estimating a moment acting on at least one joint of each leg body of the bipedal walking moving body 
by the use of the estimated values of the floor reaction force on each leg body sequentially derived by the method 
according to claim 4, comprising: 

a step of sequentially deriving inclination angles, in said absolute coordinate system, of the respective rigid- 
body corresponding portions of the bipedal walking moving body corresponding to the respective rigid bodies 
of said rigid body link model, based on the inclination angle of said body, the bending angles respectively of 
said hip joint and said knee joint of each leg body, and said rigid body link model; 

a step of sequentially deriving positions of the centers of gravity of the respective rigid-body corresponding 
portions relative to said reference point, based on the inclination angles of said respective rigid-body corre- 
sponding portions, the prederived weights of the respective rigid-body corresponding portions, and the posi- 
tions of the centers of gravity of the rigid-body corresponding portions in the respective rigid-body correspond- 
ing portions, and sequentially deriving accelerations of the centers of gravity of the respective rigid-body cor- 
responding portions relative to said reference point based on time series data about the positions of the centers 



23 



1415770A1 I > 



EP 1 415 770 A1 



of gravity of the respective rigid-body corresponding portions; 

a step of sequentially deriving accelerations of the centers of gravity of the respective rigid-body corresponding 
portions in said absolute coordinate system from the accelerations of the centers of gravity of the respective 
rigid-body corresponding portions relative to said reference point, and the accelerations of said reference point 
in said absolute coordinate system; 

a step of sequentially deriving angular velocities of the respective rigid-body corresponding portions based on 
time series data about the inclination angles of the respective rigid-body corresponding portions; 
a step of sequentially deriving estimated positions of a floor reaction force acting point of each leg body in 
said bipedal walking moving body, based on at least one of an inclination angle of a thigh portion of said leg 
body and the bending angle of the knee joint of said teg body as the rigid-body corresponding portions of said 
bipedal walking moving body; and 

a step of estimating a moment acting on at least one of the joints of each leg body of said bipedal walking 
moving body based on an inverse dynamics model using the estimated value of said floor reaction force, the 
estimated position of said floor reaction force acting point, the accelerations of the centers of gravity of the 
respective rigid-body corresponding portions and the angular velocities of said rigid-body corresponding por- 
tions in said absolute coordinate system, the inclination angles of the respective rigid-body corresponding 
portions, the prederived weights and sizes of the respective rigid-body corresponding portions, the positions 
of the prederived centers of gravity of the rigid-body corresponding portions in the respective rigid-body cor- 
responding portions, and prederived moments of inertia of the respective rigid-body corresponding portions. 
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